Systems and Methods of Tracking Locations Visited by Mobile Devices to Quantify a Change Computed based on Matching Populations used in Change Measurement

ABSTRACT

Systems and methods including mobile devices determining their locations using location determination units, such as GPS receivers. A computing device generates a statistical measurement of location pattern changes in relation to a predetermined region for a predetermined period of time. The measurement is generated from a difference between responses of exposed mobile devices and responses of non-exposed mobile devices, which devices are selected based on a cluster analysis to match the distributions of the exposed mobile devices and the non-exposed mobile devices in an attribute space. Each response of a mobile device can be computed based on applying, at a time instance, an anti-symmetric weight function having an exponential distribution to a time series of visitation measurements to generate a response corresponding to the time instance, and summing the responses corresponding to time instances falling within the predetermined period of time as the response of the mobile device.

RELATED APPLICATIONS

The present application is a continuation application of U.S. patent application Ser. No. 15/990,521 filed May 25, 2018, issued as U.S. Pat. No. 10,779,109 on Sep. 15, 2020, and entitled “Systems and Methods of Tracking Locations Visited by Mobile Devices to Quantify a Change Computed based on Matching Populations used in Change Measurement”, which is a continuation application of U.S. patent application Ser. No. 15/615,669 filed Jun. 6, 2017, issued as U.S. Pat. No. 10,323,944 on Jun. 18, 2019, and entitled “Systems and Methods of Tracking Locations Visited by Mobile Devices to Quantify a Change Computed based on Matching Populations used in Change Measurement”, which claims the benefit of the filing date of Prov. U.S. Pat. App. Ser. No. 62/346,689, filed Jun. 7, 2016, where the present application is also a continuation-in-part application of U.S. patent application Ser. No. 15/474,915 filed Mar. 30, 2017, issued as U.S. Pat. No. 10,036,638 on Jul. 31, 2018, and entitled “Systems and Methods of Tracking Locations Visited by Mobile Devices Residing within a Predetermined Activity Area to Quantify a Change”, which claims the benefit of the filing date of Prov. U.S. Pat. App. Ser. No. 62/346,689, filed Jun. 7, 2016; U.S. patent application Ser. No. 15/435,287 filed Feb. 16, 2017, issued as U.S. Pat. No. 9,986,527 on May 29, 2018, and entitled “Systems and Methods of Tracking Locations Visited by Mobile Devices to Quantify a Change from a Time Series of Responses”, which claims the benefit of the filing date of Prov. U.S. Pat. App. Ser. No. 62/429,274 filed Dec. 2, 2016 and Prov. U.S. Pat. App. Ser. No. 62/346,689 filed Jun. 7, 2016; U.S. patent application Ser. No. 15/435,185 filed Feb. 16, 2017, issued as U.S. Pat. No. 10,327,094 on Jun. 18, 2019, and entitled “Systems and Methods to Track Locations Visited by Mobile Devices and Determine Neighbors of and Distances among Locations”, which claims the benefit of the filing dates of Prov. U.S. Pat. App. Ser. No. 62/377,256 filed Aug. 19, 2016 and Prov. U.S. Pat. App. Ser. No. 62/346,689 filed Jun. 7, 2016; U.S. patent application Ser. No. 15/697,342 filed Sep. 6, 2017, issued as U.S. Pat. No. 10,142,788 on Nov. 27, 2018, and entitled “Systems and Methods to Track Regions Visited by Mobile Devices and Detect Changes in Location Patterns”, which is a continuation application of U.S. patent application Ser. No. 15/174,608 filed Jun. 6, 2016, and issued as U.S. Pat. No. 9,769,619 on Sep. 19, 2017, which is a continuation application of U.S. patent application Ser. No. 14/679,909 filed Apr. 6, 2015, and issued as U.S. Pat. No. 9,374,671 on Jun. 21, 2016; U.S. patent application Ser. No. 15/655,751 filed Jul. 20, 2017, issued as U.S. Pat. No. 10,034,141 on Jul. 24, 2018, and entitled “Systems and Methods to Identify Home Addresses of Mobile Devices”, which is a continuation application of U.S. patent application Ser. No. 15/337,939 filed Oct. 28, 2016, and issued as U.S. Pat. No. 9,730,028 on Aug. 8, 2017, which is a continuation application of U.S. patent application Ser. No. 15/014,943 filed Feb. 3, 2016, and issued as U.S. Pat. No. 9,503,863 on Nov. 22, 2016, which is a continuation application of U.S. patent application Ser. No. 14/594,010 filed Jan. 9, 2015, and issued as U.S. Pat. No. 9,291,700 on Mar. 22, 2016, the entire disclosures of which applications are hereby incorporated herein by reference.

The present application relates to U.S. patent application Ser. No. 15/435,247, filed Feb. 16, 2017 and entitled “Systems and Methods to Track Locations Visited by Mobile Devices and Update Keyword Values based on Neighbor Relations and Distances among Locations”, and U.S. patent application Ser. No. 14/593,947, filed Jan. 9, 2015 and issued as U.S. Pat. No. 9,307,360 on Apr. 5, 2016, which has a continuation application Ser. No. 15/016,067, filed Feb. 4, 2016 and published as U.S. Pat. App. Pub. No. 2016/0205503. The present application also relates to U.S. Pat. App. Pub. Nos. 2014/0012806, and 2015/0052132. The entire disclosures of the above identified patents and/or patent applications are hereby incorporated herein by reference.

FIELD OF THE TECHNOLOGY

At least one embodiment of the disclosure relates to computational efficient ways to identify location visitation measurements and measure a change in location behaviors of mobile devices.

BACKGROUND

A location determination system, such as a Global Positioning System (GPS), allows a mobile device, such as a mobile phone, a smart phone, a personal media player, a GPS receiver, etc., to determine its current location on the earth. The location of the mobile device is typically calculated as a set of coordinates, such as the longitude and latitude coordinates of a point on the surface of the earth.

However, the location of the mobile device in the form of coordinates of a point on the surface of the earth does not provide sufficient information of interest about the location, such as whether the mobile device is within a particular region associated with a set of known properties.

For example, it may be of interest in certain applications to determine whether the location of the mobile device is within the store of a merchant, within the home of the user of the mobile device, within a recreation area, within a commercial district, etc.

For example, U.S. Pat. App. Pub. No. 2014/0012806, published Jan. 9, 2014 and entitled “Location Graph Based Derivation of Attributes”, discusses the generation of a user profile based on mapping the locations of a mobile device to predefined geographical regions and use the attributes associated with the predefined geographical regions to derive and/or augment the attributes of the user profile.

For example, U.S. Pat. App. Pub. No. 2008/0248815, published Oct. 9, 2008 and entitled “Systems and Methods to Target Predictive Location based Content and Track Conversions”, discusses the need to analyze the location of a mobile device to determine the types of businesses that the user of the mobile device typically visits, or visited. When the location of a mobile device is within a predefined distance from either the address of a particular business or a geographic location associated with the business, or within a geometric perimeter of the particular business location, it may be determined that the user of the mobile device was at the particular business.

Ray Casting is a known technology to determine whether a given point is within a polygon represented by a set of vertexes. However, Ray Casting is computational intensive involving floating point number computations.

The Military Grid Reference System (MGRS) is a standard used for locating points on the earth. It uses grid squares of various lengths at different resolutions, such as 10 km, 1 km, 100 m, 10 m, or 1 m, depending on the precision of the coordinates provided. An MGRS coordinate includes a numerical location within a 100,000 meter square, specified as n+n digits, where the first n digits give the easting in meters, and the second n digits give the northing in meters.

The disclosures of the above discussed patent documents are hereby incorporated herein by reference.

BRIEF DESCRIPTION OF THE DRAWINGS

The embodiments are illustrated by way of example and not limitation in the figures of the accompanying drawings in which like references indicate similar elements.

FIG. 1 shows a system to determine whether a mobile device is within a region having a predetermined geographical boundary according to one embodiment.

FIGS. 2-4 illustrate a grid system used to determine whether a location of a mobile device is within the geographical boundary of a region according to one embodiment.

FIGS. 5-7 illustrate a hierarchical grid system used to determine whether a location of a mobile device is within the geographical boundary of a region according to one embodiment.

FIGS. 8 and 9 show a top level grid and the identification of cells within the grid according to one embodiment.

FIG. 10 shows an intermediate level grid and the identification of cells within the grid according to one embodiment.

FIG. 11 shows the identification of cells within a grid having the finest resolution in a grid hierarchy according to one embodiment.

FIG. 12 shows the method to determine whether a location of a mobile device is within the geographical boundary of a region according to one embodiment.

FIG. 13 illustrates an example of converting the coordinates of a location to an identifier of a cell and converting the identifier of the cell to the coordinates of a vertex of the cell according to one embodiment.

FIG. 14 shows a system configured to map a location of a mobile device to one or more identifications of regions according to one embodiment.

FIG. 15 illustrates a data processing system according to one embodiment.

FIG. 16 shows a method of mapping a location of a mobile device to a region according to one embodiment.

FIG. 17 shows a method to detect differences in location patterns according to one embodiment.

FIG. 18 shows a method to detect differences in location patterns of different mobile devices visiting a predetermined region according to one embodiment.

FIG. 19 shows a method to measure the influence of an event based on differences in location patterns of mobile devices visiting a predetermined region according to one embodiment.

FIG. 20 shows a method to identify mobile devices having similar patterns of visiting a predetermined region according to one embodiment.

FIG. 21 shows a method to identify mobile devices having similar patterns of visiting predetermined regions according to one embodiment.

FIG. 22 shows a method to measure the influence of an event according to one embodiment.

FIGS. 23-25 illustrate a system to organize location data via a grid system according to one embodiment.

FIG. 26 illustrates a location data processing system to establish a graph of locations according to one embodiment.

FIG. 27 shows a method to generate a location graph according to one embodiment.

FIGS. 28 and 29 illustrate a method to propagate keyword values from neighboring locations based on distances to neighboring locations according to one embodiment.

FIGS. 30 and 31 illustrate a method to propagate keyword values from neighboring locations based on visitation counts of locations according to one embodiment.

FIGS. 32 and 33 illustrate a method to propagate keyword values from mobile devices to a location visited by the mobile devices according to one embodiment.

FIGS. 34 and 35 illustrate a method to propagate keyword values from locations visited by a mobile device to the mobile device according to one embodiment.

FIG. 36 shows a method to update keyword values according to one embodiment.

FIG. 37 illustrates a method to construct a visitation measurement for a time interval according to one embodiment.

FIG. 38 illustrates a method to construct a location response at a time instance from visitation measurements of a plurality of time interval according to one embodiment.

FIG. 39 illustrates a method to construct a location response for a time period according to one embodiment.

FIG. 40 shows a method to determine a measurement of a location change according to one embodiment.

FIG. 41 shows a method to measure a statistical change in location behavior of mobile devices according to one embodiment.

FIG. 42 illustrates an activity area for the selection of mobile devices for measuring a statistical change in location behavior of mobile devices according to one embodiment.

FIG. 43 shows a method to measure a statistical change in location behavior of mobile devices based on an activity area according to one embodiment.

FIG. 44 shows a method to identify an activity area.

FIG. 45 shows another method to identify an activity area.

FIG. 46 shows a further method to identify an activity area.

FIG. 47 illustrates clustering mobile devices in an attribute space to match exposed mobile devices and non-exposed mobile devices.

FIG. 48 illustrates a selection of unbalanced matching of exposed and non-exposed mobile devices in clusters.

FIG. 49 illustrates a selection of balanced matching of exposed and non-exposed mobile devices in clusters.

FIG. 50 shows a method to compute a visitation change measurement using balanced matching of mobile devices.

FIG. 51 shows a method to compute a visitation change measurement using balanced or unbalanced matching of mobile devices.

FIG. 52 shows a method to compute a visitation change measurement based on matching population distributions.

DETAILED DESCRIPTION

The following description and drawings are illustrative and are not to be construed as limiting. Numerous specific details are described to provide a thorough understanding. However, in certain instances, well known or conventional details are not described in order to avoid obscuring the description. References to one or an embodiment in the present disclosure are not necessarily references to the same embodiment; and, such references mean at least one.

One embodiment of the disclosure provides a computationally efficient method and system to determine whether a location of the mobile device is within a predetermined geographical boundary of a region and/or to determine, among a plurality of predefined regions, the identity of one or more regions within which the location of the mobile device is positioned.

FIG. 1 shows a system to determine whether a mobile device is within a region having a predetermined geographical boundary according to one embodiment.

In FIG. 1, a location determination system uses the wireless signals (e.g., 179) transmitted to and/or from the mobile device (109) to determine the location (111) of the mobile device (109) on the surface of the earth.

For example, the location determination system may use Global Positioning System (GPS) satellites (e.g., 117) (and/or base stations (e.g., 115)) to provide GPS signals to the mobile device (109). The mobile device (109) is configured to determine the location (111) of the mobile device (109) based on the received GPS signals. In general, multiple GPS satellites (e.g., 117) and/or base stations (e.g., 115) are used to provide the wireless signals (e.g., 179) from different locations for a GPS receiver to determine its locations.

In FIG. 1, the mobile device (109) is configured with a cellular communications transceiver to communicate with the base stations (e.g., 113, 115) of a cellular communications network.

For example, in one embodiment, the mobile device (109) is configured to use signal delays in the cellular communications signals to or from a plurality of cellular base stations (e.g., 113, . . . , 115) to compute the location coordinates of the mobile device (109).

In FIG. 1, a server (187) is configured to communicate with the mobile device (109) via the network (189) and the cellular communications infrastructure (e.g., the base station (113)). The server (187) is connected to a database (181) storing information about the predefined regions (e.g., 101, 103, . . . 105, 107).

For example, the database (181) is configured to store the identifications of a set of cells that are within the boundary of a region (e.g., 101). The server (187) is configured to convert the location (111) of the mobile device (109) to a cell identification and search the identifications of the set of cells representing the region (101) to determine if the cell identification converted from the location (111) of the mobile device (109) is in the set of cell identifications representing the region (101). If the cell identification of the location (111) is found in the set of cell identifications representing the region (101), the location (111) is considered being within the boundary of the region (e.g., 101).

In one embodiment, a hierarchical grid system is used to construct cells that are within the boundary of the region (e.g., 101). Thus, the number of cells within the region (e.g., 101) can be reduced, while the precision of the determination can be selected at a desired level (e.g., 1 meter).

In one embodiment, the identifications of the cells are configured to be signed integer numbers. Thus, any known technologies for searching a given number within a set of signed integer numbers can be used to efficiently determine whether the cell identifier of a location (111) is within the set of cell identifiers of the region (101).

In one embodiment, the conversion of the location coordinates to a cell identifier is configured for improved computation efficiency. The cell identifier is also configured for efficient determination of the resolution of the grid in which the cell is located, the coordinates of the vertexes of the cell, the bounding boxes of the cell, and the identifications of the neighbors of the cells. Details and examples are provided below.

In one embodiment, a given region (e.g., polygon) on earth is represented by a set of cells in a hierarchical, regular grid in a longitude latitude space. In the longitude latitude space, the cells are uniform rectangles/squares at a given resolution; the cell identifies are constructed from the digits of the longitude/latitude coordinates for improved efficiency in conversion between coordinates and cell identifiers. In one embodiment, the resolution levels of the grids correspond to the precision of the longitude/latitude coordinates in terms of the number of digits used to after the decimal point to represent the longitude/latitude coordinates.

At a given resolution in the grid, the identity of the cell that contains a given point identified by a longitude/latitude pair can be computed via simple manipulations of the digits of the longitude/latitude pair, as illustrated in FIG. 13.

FIGS. 2-4 illustrate a grid system used to determine whether a location of a mobile device is within the geographical boundary of a region according to one embodiment.

In FIG. 2, a grid (121) of cells is used to identify an approximation of the region (101) at a given level of resolution of the grid (121). The resolution level corresponds to the size of the cells in the grid (121).

In FIG. 2, the region (101) is represented as a polygon having a set of vertexes (e.g., 123). The set of line segments connecting the neighboring vertexes (e.g., 123) of the region (101) defines the boundary of the region (101).

FIG. 3 illustrates the selection of a set of cells (e.g., 127) that are considered to be within the boundary of the region (101). Various different methods and/or criteria can be used to classify whether a cell is within the boundary of the region (101), especially the cells that are partially in the region (101) and contain a portion of the boundary of the region (101). The disclosure of the present application is not limited to a particular way to identify or classify whether a cell that is within the boundary of the region (101).

For example, a cell may be classified as being with the region (101) when the overlapping common portion between the cell and the region (101) is above a predetermined percentage of the area of the cell.

For example, a cell may be classified as being with the region (101) when a length of one or more segments of the region (101) going through the cell is above a threshold.

For example, the vertexes of the region (101) may be mapped to the nearest grid points to determine an approximation of the boundary of the region (101) that aligns with the grid lines to select the cells that are located within the approximated boundary of the region (101).

FIG. 4 illustrates the determination of the location (111) within the set of cells (131, . . . , 133, . . . , 139) according to one embodiment. In FIG. 4, each of the cells (131, . . . , 133, . . . , 139) represents a portion of the region (101). To determine whether the location (111) is within the boundary of the region (101), the system is configured to determine whether the set of cells (131, . . . , 133, . . . , 139) contains the location (111).

In one embodiment, to efficiently determine whether any of the cells (131, . . . , 133, . . . , 139) contains the location (111), each of the cells (131, . . . , 133, . . . , 139) is assigned a cell identifier. In one embodiment, each of the cell identifier is a signed integer for improved computation efficiency; and the cell identifier is configured in such a way that the coordinates of any location within the cell can be manipulated via a set of predetermined, computationally efficient rules to provide the same cell identifier, as further illustrated in FIGS. 12 and 13.

In FIG. 4, after the coordinates of the location (111) is converted to the cell identifier of the cell (133) that contains the location (111), the system determines whether the location (111) is within the region corresponding to the set of cells (131, . . . , 133, . . . , 139) by searching in the cell identifiers of the set of cells (131, . . . , 133, . . . , 139) representative of the region (101) to find a match to the cell identifier of the cell (133) that is converted from the coordinates of the location (111). If a match is found, the location (111) is determined to be within the region (101); otherwise, the location (111) is determined to be outside of the region (101).

To improve the accuracy in the approximation of the region (101) and computational efficiency, the cells of a hierarchical grid system is used to approximate the region (101). FIGS. 5-7 illustrate a hierarchical grid system used to determine whether a location of a mobile device is within the geographical boundary of a region according to one embodiment.

In FIG. 5, grids of different resolutions are used to identify a set of cells to approximate the region (101). The grids has a predetermined hierarchy, in which the grid lines of a high level grid aligns with some of the grid lines of a low level grid such that the cells of the low level grid subdivide the cells of the high level grid. The grids of different resolutions have different cell sizes.

In general, a grid having a higher resolution and thus smaller cell size can approximate the region (101) in better precision, but uses more cells.

In one embodiment, the cells from the lower resolution grid is used in the interior of the region (101) to reduce the number of cells used; and the cells from the higher resolution grid is used near the boundary of the region (101) to improve precision in using the set of cells to approximately represent the region (101).

For example, in one embodiment, the lowest resolution gird is applied to identify a set of cells to approximate the region (101). The cells in the lowest resolution grid that contain the boundary of the region (101) are split in accordance with the grid of the next resolution level to identify cells in the grid of the next resolution level for improved precision in representing the region (101). The cell splitting process can be repeated for further improved precision using a higher resolution grid.

FIG. 6 illustrates the use of cells from two levels of hierarchical grids to approximate the region (101).

After the set of cells used to approximate the region (101) are identified (e.g., as illustrated FIG. 6), the system is configured to determine whether the location (111) of the mobile device (109) is within the region (101) based on whether any of the set of cells representing the region contains the location (111), in a way as illustrated in FIG. 7.

For example, in one embodiment, each of the cells used in FIG. 7 to represent a part of the region (101) is provided with a cell identifier. The coordinates of the location (111) is mapped to a cell identifier at a given resolution level. The system is configured to search in the set of cell identifiers of region (101) at the corresponding resolution level to determine whether there is a match to the cell identifier as determined from the coordinates of the location (111). If a match in cell identifier is found at any resolution level, the location (111) is determined to be within the region (101) represented by the set of cells; otherwise, the location (111) is determined to be outside the boundary of the region (101).

In one embodiment of FIG. 1, a hierarchical grid system is used to approximate the predefined regions (101, 103, . . . , 105, 106) with cells. Each of the cells is classified/identified as being in one or more of the regions (101, 103, . . . , 105, 106). The database (181) stores the identifiers of the cells in association with the identifies of the respective regions (101, 103, . . . , 105, 106); and the server (187) is configured to compute the identifiers of the cells of different resolutions that contain the location (111) and determine if any of the cell identifiers stored in the database (181) in association with the identifiers of the regions (101, 103, . . . , 105, 106) has the same cell identifier as the location (111). If a matching cell identifier is found, the location (111) of the mobile device (109) is determined to be with the respective region(s) (e.g., 101) associated with the corresponding cell identifier; otherwise, the location (111) is determined to be outside all of the regions (101, 103, . . . , 105, 106) represented by the set of cell identifiers stored in the database (181).

FIGS. 8 and 9 show a top level grid and the identification of cells within the grid according to one embodiment.

In one embodiment, the location (111) of the mobile device (109) is determined to be on the surface of the earth in terms of the longitude and latitude coordinates. In a coordinate system as illustrated in FIG. 8, the longitude coordinates are configured to be within the range of −180 degrees to 180 degrees; and the latitude coordinates are configured to be with the range of −90 degrees to 90 degrees.

In one embodiment, a hierarchical grid system on the surface of the earth is based on a regular grid in the longitude latitude space illustrated in FIG. 9.

In FIG. 9, the cells in the top level grid have a uniform size of 10 degrees in longitude and 10 degrees in latitude. In FIG. 9, the cells are identified by the row identifiers ranging from −9 to −1 and 1 to 9 and column identifiers ranging from 1 to 36.

In FIG. 9, the row and column identifiers are configured in a way to avoid using zero as a row identifier or a column identifier.

In FIG. 9, the row identifier of 1 is assigned to the row of cells between 0 to 10 degrees of latitude; the row identifier of 2 is assigned to the row of cells between 10 to 20 degrees of latitude; etc. The rows of cells between 0 to −90 degrees of latitudes are assigned similar row identifiers with a negative sign. For example, the row identifier of −1 is assigned to the row of cells between 0 to −10 degrees of latitude; the row identifier of −2 is assigned to the row of cells between −10 to −20 degrees of latitude; etc. As a result, the row identifier has a sign and a single digit for the top level cells illustrated in FIG. 9; and the single digit is not zero for any of the rows. Thus, for each location that is inside a cell in the top level grid as illustrated in FIG. 9, the row identifier of the cell containing the location has the same sign as the latitude coordinate of the location and the single digit that equals to 1 plus the tens digit of the latitude coordinate of the location.

In FIG. 9, the column identifier of 1 is assigned to the column of cells having longitude coordinates between −180 to −170 degrees; the column identifier of 2 is assigned to the column of cells having longitude coordinates between −170 to −160 degrees; etc. Thus, for each location that is inside a cell in the top level grid as illustrated in FIG. 9, the column identifier of the cell containing the location has no sign (e.g., the column identifier is always greater than zero) and one or two digits that correspond to adding 18 to a number formed by using the hundreds digit of the longitude as the tens digit and the tens digit of the longitude as the ones digit.

The combination of the row identifier and the column identifier of a cell uniquely identifies the cell within the top level grid as illustrated in FIG. 9. For example, the digits of the column identifier can be appended to the row identifier to generate a signed number that uniquely identifies the cell within the grid illustrated in FIG. 9. For a given cell identifier, the row identifier and the column identifier can be unambiguously deduced from the cell identifier itself, since the row identifier has a signal digit and a sign. The longitude and latitude coordinates of the vertexes of the cell can be computed from the row identifier and the column identifier.

Although FIG. 9 illustrates a preferred way to code the row identifiers and the column identifiers based on the longitude and latitude coordinates of the locations within the cells, alternative coding schemes can be used.

For example, the rows can be coded from 1 to 18 for latitudes from −90 degrees to 90 degrees; and the columns can be coded from 10 to 45 for longitudes from −180 degrees to 180 degrees. Thus, both the row and column identifiers are positive integers, while the column identifiers always have two digits.

For example, the rows can be coded from 11 to 28 for latitudes from −90 degrees to 90 degrees; and the columns can be coded from 11 to 46 for longitudes from −180 degrees to 180 degrees. Thus, both the row and column identifiers are positive integers having two digits.

FIG. 10 shows an intermediate level grid and the identification of cells within the grid according to one embodiment. In FIG. 10, a given cell at a higher level grid (e.g., a cell in the top level grid as illustrated in FIG. 9) is subdivided into 10 rows and 10 columns. The coding of the rows and columns correspond to the measurement directions of the longitude and latitudes coordinates such that the corresponding digits in the longitude and latitudes coordinates at a given precision level can be used directly as the row and column identifiers of the sub-cells within the cell at the higher level grid.

For example, when the cell that is being subdivided into the 10 rows and 10 columns has a size of 10 degrees in longitude and 10 degrees in latitude (e.g., as illustrated in FIG. 9), the row identifier and column identifier of the sub-cells correspond to the ones digit of the latitude and longitude coordinates of the points within the respective sub-cells.

For example, when the cell that is being subdivided into the 10 rows and 10 columns has a size of 1 degree in longitude and 1 degree in latitude, the row identifier and column identifier of the sub-cells correspond to the one-tens digit of the latitude and longitude coordinates of the points within the respective sub-cells.

FIG. 11 shows the identification of cells within a grid having the finest resolution in a grid hierarchy according to one embodiment. In FIG. 11, the row identifiers and column identifiers are padded by 1, in comparison with the row and column coding scheme illustrated in FIG. 10.

In one embodiment, an identifier cell for a given resolution includes sufficient information to identify the corresponding cells in the higher level grid(s) that contains the cell. Thus, a cell identifier uniquely identifies a cell in the entire hierarchical grid without ambiguity.

FIG. 12 shows the method to determine whether a location of a mobile device is within the geographical boundary of a region according to one embodiment.

In FIG. 12, the location (111) of the mobile device (109) is determined in terms of the longitude coordinate (143) and the latitude coordinate (145).

For a given resolution level (147), the longitude coordinate (143) and the latitude coordinate (145) are truncated to generate the column identifier (149) and the row identifier (151). Applying (155) the resolution level (147) includes truncating the longitude coordinate (143) and the latitude coordinate (145) to the corresponding digits of precision to generate the column identifier (149) and the row identifier (151). In one embodiment, the digits corresponding to the top level grid and the bottom level grid at the given resolution are adjusted according to FIGS. 9 and 11.

In FIG. 12, the column identifier (149) and the row identifier (151) are combined to generate the cell identifier (153) of the location (111) of the mobile device at the given resolution level (147).

In one embodiment, the database (181) stores a set of cell identifiers (161, . . . , 163) that are associated with the region (101) defined by a predetermined boundary. The server (187) searches (157) the set of cell identifiers (161, . . . , 163) to find a match with the cell identifier (153). If a match is found, the location (111) of the mobile device (109) is determined to be within the boundary of the region (101).

In one embodiment, the database (181) stores a set of cell identifiers (e.g., 161, . . . , 163, 165, . . . ) associated with respective different regions (e.g., 101, 103, . . . ). When the cell identifier (153) of the location (111) of the mobile device (109) is found to be matching with a particular cell identifier (e.g., 163 or 165), the region (e.g., 101 or 103) associated with the particular cell identifier (e.g., 163 or 165) is determined to be the region in which the mobile device (141) is located.

In one embodiment, when a cell contains the boundary of two regions (e.g., 101 and 103), the cell identifier of the cell can be associated with both regions (e.g., 101 and 103). The system may optionally further determine which region the cell is in based on the coordinates of the vertexes defining the boundary (or other parameters that define the boundary between the regions).

FIG. 13 illustrates an example of converting the coordinates of a location to an identifier of a cell and converting the identifier of the cell to the coordinates of a vertex of the cell according to one embodiment.

In FIG. 13, the location has a latitude coordinate of −51.12345678 and a longitude coordinate of −41.12345678. A resolution at the fifth digit after the decimal point is applied to the coordinates to generate the truncated coordinates (−41.12345, −51.12345). The decimal point is removed to obtain the longitude digits −4112345 and the latitude digits −5112345. Since the length of the equator of the earth is about 40,075 km, the cell size near the equator is about 1.11 meters at the resolution corresponding to the fifth digit.

In accordance with the scheme for the top level grid illustrated in FIG. 9, the tens digit for the latitude coordinate is padded with one (without considering the sign of the latitude); and the hundreds digit and tens digit, including the sign, of the longitude coordinate is padded with 18 to generate the row identifier −6 and the column identifier 14 for the top level grid.

In accordance with FIG. 10, the row identifiers and column identifiers of the sub-cells in the hierarchical grid correspond to the respective latitude digits and longitude digits (1, 1, 2, 3, 4).

In accordance with FIG. 11, the row identifiers and column identifiers of the sub-cells in the bottom hierarchy is padded with 1, if the longitude and/or the latitude coordinates of the location is not exactly on the grid lines of the resolution level (e.g., if the longitude or latitude coordinate has nonzero digits after the fifth digit behind the decimal point). One is not padded at the last digit when the longitude and/or the latitude coordinates of the location is exactly on the grid lines of the resolution level (e.g., if the longitude or latitude coordinate has no nonzero digits after the fifth digit behind the decimal point). According to this padding scheme, in the northern hemisphere locations on the northern edge of a cell are included in the cell but not the locations on the southern edge. In the southern hemisphere, locations on the southern edge of a cell are included in the cell but not the locations on the northern edges. Locations on the eastern edge of a cell are included in the cell, but not the western edge.

Thus, the location (−41.12345678, −51.12345678) has the row and column identifiers −6112346 and 14112346. The digits of the column identifier are appended to the digits of the row identifier to generate the cell identifier −611234614112346.

In FIG. 13, the row and column identifiers can be recovered from the cell identifier. The number of digits in the cell identifier divided by 2 provides the number of leading digits for the row identifier; and the remaining digits are for the column identifier. From the row identifier and column identifiers, the latitude digits and longitude digits can be computed via subtraction of the respective padding. The truncated coordinates can be computed from the latitude digits and longitude digits respectively, which can be used to determine the coordinates of a vertex of the cell as (−41.12345, −51.12345). Based on the resolution of the cell being at 0.00001, the coordinates of other vertexes of the cell can be determined as (−41.12346, −51.12345), (−41.12346, −51.12344), (−41.12345, −51.12344). The bounding box of the cell and the neighboring cells can also be easily identified based on the coordinates.

FIG. 13 illustrates a way to append the digits of the column identifier to the digits of the row identifier to generate the cell identifier. Alternatively, the row identifier and the column identifier can be combined in other ways that can be reversed to derive the row identifier and the column identifier from the cell identifier.

For example, when the top level column identifiers are mapped to the range 11 to 46 to have a fixed number of two digits for the top level column, the column identifier is 2411236. Since there is no ambiguity in the number of digits used to represent the top level column, the top level column identifier (24) can be appended after the top level row identifier (−6), which is then appended with the row and column identifiers of the next level, and so on. Thus, a cell identifier of −6241111223366 can be generated, with the sign then the first three digits representing the top level row and column, and two digits for subsequent next level row and column to identifying the subdivision within the higher level cell.

In some embodiments, the row and column identifiers of the bottom level are not padded in a way illustrated in FIG. 11 to have different ways to account for the locations on grid lines at the lowest level resolution.

FIGS. 9-11 and 13 illustrate a grid hierarchy based on a decimal representation of longitude and latitude coordinates. Alternatively, the grid hierarchy can be constructed in accordance with longitude and latitude coordinates expressed using other bases, such as binary, ternary, quintal, octal, duodecimal, etc. in a similar way.

Further, in some embodiments, the longitude and latitude coordinates may be normalized (e.g., in the standardized data range between 0 to 1); and the grids can be constructed in the space of the normalized longitude and latitude coordinates.

The hierarchical grid can also be extended to a three-dimensional space. For example, a hierarchical grid can be constructed with regular grids in the longitude, latitude, altitude space, or in a mapped or normalized longitude, latitude, and altitude space.

FIG. 14 shows a system configured to map a location of a mobile device to one or more identifications of regions according to one embodiment. In FIG. 14, the mobile device (109) determines the coordinates (171) of its location (111) based on the wireless signals (179) to and/or from a location determination system, such as the Global Positioning System (GPS).

The coordinates (171) are converted to a cell identifier (173) of a cell that contains the location, e.g., in a way as illustrated in FIG. 12 or 13.

In the database (181), a set of cell identifiers are stored in association with region identifiers (185), where each of the cell identifiers is associated with one or more of the respective regions when the respective cell contains at least a portion of the one or more of the respective regions.

In one embodiment, the set of cell identifiers are organized as a cell identifier tree (183) to facilitate the search of a matching identifier.

For example, the cell identifier tree (183) can be constructed as a self-balancing tree for efficient searching of a cell identifier matching the cell identifier (173) generated from the coordinates (171) of the mobile device (109).

In general, any methods to search for an identifier with a set of predetermined identifiers can be used to search for the matching cell identifier (173).

From the association of the cells with the region identifiers (185) in the database, the server (187) determines the identification (175) of the one or more defined regions that are at least partially in the cell identified by the cell identifier (173). Thus, the location (111) of the mobile device (109) is determined to be within the region(s) identified by the identification (175) of the defined region(s).

Similarly, after regions of different sizes and locations are represented via the cells in the hierarchical grid, the system can be configured to efficiently compute overlapping portions of regions via searching for cells having the same identifications.

For example, to determine the approximate overlapping between regions, the percentage of overlapping, the square of overlap, etc., the system is configured to count a number of overlapped cells to determine the overlapping.

In one embodiment, a polygon or any other shape is approximated by a set of rectangular and/or square cell of different sizes in a suitable coordinate system (e.g., in longitude latitude space). Each cell is represented by a single number as identifier. The identifiers of the cells used to approximate the polygon or shape can be organized as a binary tree, a self-balanced tree, a Red/Black Tree, or other structures that are known to provide logarithmic search time to improve the computation efficiency in determining whether a point is within the polygon or shape.

For example, a polygon representing the boundary of United States of America USA on a map may include 2,000 vertexes. The Ray Casting algorithm has O(n) complexity to calculate if a point is within the polygon. When this polygon is approximated via a hierarchical grid system discussed above, the polygon can be represented 700 to 2,000,000 cells in the longitude latitude space, depending on the required precision. When the polygon is represented by 2,000,000 cells and their corresponding identification numbers, searching a matching identifier at the same precision via a binary tree gives log(2,000,000)=21 complexity, which is much less than 2,000. Thus, the present disclosure improves the computational efficiency of identifying a region in which a mobile device is located.

FIG. 16 shows a method of mapping a location of a mobile device to a region according to one embodiment. For example, the method of FIG. 16 can be implemented in the system of FIG. 1 and/or FIG. 14, using the grid system illustrated FIGS. 2-8, and/or the grid system and cell identifier system illustrated in FIGS. 8-13.

In FIG. 16, a computing apparatus is configured to: identify (221) a set of cells in a grid system that are within the predefined boundary of a geographic region; receive (223) a location (111) of a mobile device (109); convert (225) the location (111) to the identifier of a cell that contains the location; and search (227) identifiers of the set of cells to determine if the cell identifier of the location (111) is in the set. If it is determined (228) that the cell identifier of the location (111) is in the set, the computing apparatus determines (229) that the location (111) of the mobile device (109) is in the geographic region.

In one embodiment, the computing apparatus includes at least one of: the database (181) and the server (187).

In one embodiment, the database (181) is configured to store an identifier of a geographical region (101) having a predefined geographical boundary defined by a set of vertexes (e.g., 123) or a set of other parameters, such as a center location and a radius.

The database (181) further stores a set of cell identifiers, each of which identifies a cell that is determined to be within the predefined geographical boundary of the geographical region (101).

After the server (187) receives, from a mobile device (109), a location (111) of the mobile device (109), the server (187) converts a set of coordinates (143, 145) of the location (111) of the mobile device (109) to a cell identifier (153) of a cell that contains the location (111). In some embodiments, the mobile device (109) generates the cell identifier (153) at a desired precision level to represent the location (111) of the mobile device (109).

The server (187) determines whether the location (111) of the mobile device (109) is within the geographical region (101) based on searching the set of cell identifiers to determine if the set has the cell identifier (153) computed from the coordinates (143, 145) of the location (111) of the mobile device (109).

In one embodiment, to convert the set of coordinates (143, 145) of the location (143, 145) to the cell identifier (153), the server (187) (or the mobile device (109)) generates two integers from longitude and latitude coordinates of the location (111) of the mobile device (109) according to a precision level (e.g., resolution level (147), and combine the two integers into the first cell identifier (153) without using a floating point number computation.

In one embodiment, each cell using the in the system to approximate the regions and the locations is a rectangle/square area in a longitude latitude space of locations on the earth. The size of the cell can be unambiguously determined from the cell identifier itself. Further, the longitude and latitude coordinates of corners of the cell identified by the cell identifier can be unambiguously determined from the cell identifier itself.

In one embodiment, the set of cells identified by the set of cell identifiers to approximate one or more regions (e.g., 101, 103, . . . , 105, . . . , 107) has a plurality of different cell sizes that correspond to a plurality of predetermined cell resolution levels. Each of the plurality of predetermined cell resolution levels corresponds to a predetermined precision level of longitudes and latitudes of locations on the earth. For example, each of the plurality of predetermined cell resolution levels corresponds to a precision to a predetermined digit after the decimal point in longitude and latitude coordinates of locations on the earth.

In one embodiment, a cell identifier itself includes sufficient information to determine the resolution level of the cell, the coordinates of the vertexes of the cell, and the identifiers of the neighboring cells, etc.

In one embodiment, the database (181) stores data mapping each cell identify in the set of cell identifiers to at least one region identifier, where the cell contains a least a part of each of the regions identified by the at least one region identifier. The server (187) is configured to search the set of cell identifiers to find a cell identifier that matches with the cell identifier (153) computed from the location (141) and thus determine at least one region identifier associated with the matching cell identifier.

For example, in one embodiment, the set of coordinates of the location (111) includes longitude (143) and latitude (145) of the location (111). To converting the coordinates (143, 145) to the cell identifier (153), the server (187) (or the mobile device (109)) selects digits from the longitude (143) and the latitude (145) of the location (111) in accordance with a cell resolution level (147) and combines the digits selected from the longitude (143) and the latitude (145) of the location (111) into an integer representing the cell identifier (153) of the location (111).

As illustrated in FIG. 13, selecting the digits from the longitude and the latitude includes: selecting digits from integer part of the longitude and a first number of digits from the longitude after the decimal point of the longitude to form an integer representation of the longitude at the cell resolution level; and selecting digits from integer part of the latitude and the same first number of digits from the latitude after the decimal point of the latitude to form an integer representation of the longitude at the cell resolution level.

In one embodiment, to generate the column identifier and row identifier of the location (111), a predetermined number (e.g., one) is added to a digit of the integer representation of the latitude that corresponds to the tens digit of the latitude; and a sign is provided to the integer representation of the latitude according to the sign of the latitude.

In one embodiment, after providing a sign to the integer representation of the longitude according to the sign of the longitude, a predetermined number (e.g., eighteen) is added to digits of the integer representation of the longitude that corresponds to the hundreds digit and tens digit of the longitude, in view of the sign provided to the integer representation of the longitude.

In one embodiment, when the latitude coordinate has a non-zero portion that is discarded during the selection of the latitude digits for the integer representation of the latitude, one is added to the ones digit of the integer representation of the latitude without considering the sign of the integer representation of the latitude. When the longitude coordinate has a non-zero portion that is discarded during the selection of the longitude digits for the integer representation, one is added to the ones digit of the integer representation of the longitude without considering the sign of the integer representation of the longitude.

In one embodiment, after the server (187) receives data representing the predefined geographical boundary of the geographical region, such as the coordinates of the vertexes of a region having a polygon shape, the coordinates of the center and the radius of a region having a circular shape, etc., the server (187) identify, in a hierarchy of cell grids, the set of cell identifiers that are determined to be within the predefined geographical boundary.

In one embodiment, when the set of cells being searched having different resolutions (cell sizes), the location (111) of the mobile device (109) is converted to a plurality of cell identifiers at the corresponding resolutions; and the server (187) is configured to search a match of any of the cell identifiers at the corresponding resolutions computed from the location (111) of the mobile device (109).

For example, the identifiers of the cells of different sizes/resolutions to represent the regions can be organized in a single tree; and the identifiers of the location (111) of the mobile device (109) of corresponding sizes/resolutions can be searched concurrently or one after another to find a match.

For example, the identifiers of the cells of different sizes/resolutions to represent the regions can be organized in separate trees according to cell sizes/resolutions; and the identifiers of the location (111) of the mobile device (109) of corresponding sizes/resolutions can be searched concurrently or one after another in the respective trees for corresponding sizes/resolutions.

In one embodiment, each grid in the hierarchy of cell grids corresponds to a rectangle/square grid in longitude latitude space of locations on the earth with a predetermined resolution level that corresponds to a precision level in a floating point decimal representation of longitude and latitude coordinates.

The server (187) and/or the database (181) can be implemented as a computer apparatus in the form of a data processing system illustrated in FIG. 15.

In one embodiment, the techniques disclosed above are used to map location histories of mobile devices into visitation histories of the mobile devices to one or more pre-defined regions. The visitation patterns of the mobile devices can be used to identify mobile devices having similar behaviors and/or different behaviors in relation with locations. For example, the differences in behaviors following an event between mobile devices that are associated with the event and mobile devices that are not associated with the event can be measured as an influence of the event, if the mobile devices have the same or similar behaviors prior to the event.

For example, the mobile devices initially having similar behavior profiles may have different behaviors after some of the mobile devices are provided with an advertisement and others are not. Thus, the difference can be detected as a measurement of the influence/effectiveness of the advertisement.

FIG. 17 shows a method to detect differences in location patterns according to one embodiment.

In FIG. 17, the server (187) is configured to identify (241) a pair of mobile devices, including a first mobile device associated with an event and a second mobile device not associated with the event, to measure the influence of the event.

For the first mobile device, the server (187) is configured to: receive (243) locations of the first mobile device; convert (245) the locations to one or more cell identifiers; search (247) for regions containing the cell identifiers; and generate (249) first data representing a location pattern based on presence of the first mobile device in the regions.

For the second mobile device, the server (187) is similarly configured to: receive (253) locations of the second mobile device; convert (255) the locations to one or more cell identifiers; search (257) for regions containing the cell identifiers; and generate (259) second data representing a location pattern based on presence of the second mobile device in the regions.

The server (187) identifies (251) a difference of the first data and the second data as a measurement of an influence of the event. Preferably, the pair of mobile devices have similar attributes and/or location patterns prior to the event. Alternatively, the change in location patterns of the first mobile device before and after the event can be compared to the change in location patterns of the second mobile device before and after the event to measure the influence of the event.

For example, the pair of mobile devices may be initially identified to have similar profiles in location patterns and/or other attributes, such as the demographic data of their users, the propensity scores of their users, etc. After the event of an advertisement being transmitted to the user of the first mobile device but not the user of the second mobile device, the difference in the location patterns represents a measurement of the influence/effectiveness of the advertisement in changing location/visitation patterns.

Further, the differences in other attributes associated with the pair of mobile devices can be determined as measurements of the influence/effectiveness of the advertisement with respect to corresponding attributes. Examples of such attributes include the enrollment in a program or service, visitation to certain websites, foot traffic to a store, web traffic to a website, etc.

For example, the event may include an advertisement for a geographical region; and the method may be adapted to measure the effectiveness of the advertisement in changing a pattern of visitation to the geographical region, as illustrated in FIG. 18.

FIG. 18 shows a method to detect differences in location patterns of different mobile devices visiting a predetermined region according to one embodiment.

In FIG. 18, the server (187) identifies (261) a pair of mobile devices, including a first mobile device associated with an event relevant to a predetermined region and a second mobile device not associated with the event, to measure an influence of an event related to the region.

For the first mobile device, the server (187) is configured to: receive (263) locations of the first mobile device; convert (265) the locations to one or more cell identifiers; determine (267) whether the locations are in the region using the techniques discussed above based on searching for matching cell identifiers; and generate (269) first data representing a presence pattern of the first mobile device in the region.

For the second mobile device, the server (187) is similarly configured to: receive (273) locations of the second mobile device; convert (275) the locations to one or more cell identifiers; determine (277) whether the locations are in the region using the techniques discussed above based on searching for matching cell identifiers; and generate (279) second data representing a presence pattern of the second mobile device in the region.

The server (187) identifies (271) a difference of the first data and the second data as a measurement of an influence of the event. Preferably, the pair of mobile devices have similar attributes and/or location patterns prior to the event. Alternatively, the change in location patterns of the first mobile device before and after the event can be compared to the change in location patterns of the second mobile device before and after the event to measure the influence of the event.

In one embodiment, the measured influence is based on the difference between two groups of mobile devices, instead of the different between two mobile devices, to account for the pattern variations in the individual mobile devices.

FIG. 19 shows a method to measure the influence of an event based on differences in location patterns of mobile devices visiting a predetermined region according to one embodiment.

In FIG. 19, the server (187) is identifies (301) two groups of mobile devices, including a first group associated with an event relevant to a predetermined region and a second group not associated with the event, to measure the influence of the event on groups of mobile devices.

The server (187) is configured to: receive (303) locations of mobile devices in the two groups; convert (305) the locations to respective cell identifiers; determine (307) whether the locations are in the region based on whether the region has the respective cell identifiers; generate (309) first data representing a presence pattern of the mobile devices in the first group in the region; and generate (311) second data representing a presence pattern of the mobile devices in the second group in the region.

The server (187) determines (313) a difference between the first data and the second data as a measurement of the influence of the event.

In one embodiment, the server (187) determines the change of presence pattern of each mobile device, the average change of presence pattern in the first group, the average change of presence pattern in the second group, and then the difference between the average changes of presence pattern in the two groups for a measurement of the influence of the event.

In some embodiments, the presence patterns of mobile devices are determined to select groups of mobile devices having similar behaviors, including the presence patterns.

FIG. 20 shows a method to identify mobile devices having similar patterns of visiting a predetermined region according to one embodiment. The mobile devices having similar behaviors can be divided into two groups to measure the influence of an event.

In FIG. 20, the server (187) is configured to: receive (321) locations of mobile devices; convert (323) the locations to respective cell identifiers; determine (325) whether the locations are in a predetermined region based on whether the region has the respective cell identifiers; identify (327) presence patterns of the mobile devices in the region; and identify (329) a subset of the mobile devices based at least in part on similarity in presence patterns.

The server (187) divides (331) the subset of mobile devices into two groups, including a first group for association with an event relevant to a predetermined region and a second group not associated with the event. The method of FIG. 20 can then be used to measure the influence of the event.

In some embodiments, an event may have influence on the visitation patterns of a mobile device in visiting multiple regions.

FIG. 21 shows a method to identify mobile devices having similar patterns of visiting predetermined regions according to one embodiment.

In FIG. 21, the server (187) is configured to: receive (341) locations of mobile devices; convert (343) the locations to respective cell identifiers; identify (345) predetermined regions in which the locations are, based on whether the regions have the respective cell identifiers; identify (347) presence patterns of the mobile devices in the regions; identify (349) a subset of the mobile devices based at least in part on similarity in the presence patterns; and divide (351) the subset of mobile devices into two groups, including a first group for association with an event relevant to a predetermined region and a second group not associated with the event.

In one embodiment, the locations of a mobile device are received as a function of time. Thus, the visitations of the mobile devices are determined as a function of time to determine presence patterns based at least in part on time.

FIG. 22 shows a method to measure the influence of an event according to one embodiment.

In FIG. 22, the server (187) is configured to: identify (361) a set of mobile devices based at least in part on similarity in associated patterns; divide (363) the set of mobile devices into two groups, including a first group for association with an event relevant to a predetermined region and a second group not associated with the event; determine (365) a difference in associated patterns between the mobile devices in the first group and the mobile devices in the second group; and generate (367) a measurement of the influence of the event based on the difference.

In one embodiment, when the set of mobile devices initially having similar but not identical associated patterns prior to the event, the server (187) is configured determine the change of associated patterns for each mobile device before and after the event. The server (187) then determines the difference in the changes of associated patterns in the first and second groups to measure the influence of the event.

For example, the above discussed techniques can be used to measure the advertising effectiveness. In various embodiments, attribute data included in a first user profile may be used to select a second user profile that is substantially similar to the first user profile. The first user profile may include an indication of exposure to advertising content data associated with a location and the second profile may not include such an indication. For example, a first user profile may be associated with a first user that has seen an advertisement for a location and the second user profile may be associated with a second user who has not seen the ad. In various embodiments, propensity score matching and/or other approaches may be used to select a second user profile. For example, a propensity score may be generated based on the attribute data in the first user profile (e.g. demographic data, behavioral data, etc.) and the propensity score may be compared to propensity scores generated for other user profiles to select a second user profile. The second user profile may, for example, be associated with a propensity score that matches (e.g., most closely matches) the propensity score associated with the first user profile.

According to some embodiments, first behavior information (e.g., a change in number/frequency of visits to a location over a period prior to and over a period after seeing an ad related to the location) may be determined based at least in part on an association between the first user profile and a location associated with the advertising content data. Second behavior information may be determined based at least in part on an association between the second user profile and the location. An advertising effectiveness value may be generated based at least in part on the first behavior information and the second behavior information.

In one embodiment, attribute data included in a first user profile may be used to select a second user profile that is substantially similar to the first user profile. In various embodiments, attribute data may include, for example, demographic data, behavioral data, data from third-party sources, psycho-graphic data, location visit frequency patterns, shopping cart spend data (e.g., including similar products and/or categories of products), and/or any other data associated with a user. A first user profile may include a user profile for a user that has been exposed to advertising content associated with a location (e.g., an advertisement to drive foot traffic to the location). In some embodiments, attribute data included in a first user profile may be compared to attribute data associated with one or more other user profiles associated with users who have not been exposed to the advertising content. And a user profile including attributes that are substantially similar to (e.g., matches) the attributes included in the first user profile may be selected. Various approaches may be used to identify (e.g., select) matching user profiles including, for example, propensity score matching, statistical matching approaches, one-to-one matching, and/or any other any other matching technique.

In various embodiments, the first user profile may include a user profile from an exposed/test group, and the second user profile may include a user profile from a control group and/or general population group. The first user profile and second user profile may be used to test (e.g., measure) the effectiveness and/or influence of advertising content data associated with a location (e.g., an advertisement to drive users to a retail location). The first user profile (e.g., the exposed/test group user profile) may include an indication that a user associated with the first user profile has been exposed to advertising content data associated with a location. And the second user profile (e.g., control group user profile, general population user profile) may include an indication that a user associated with the second user profile has not been exposed to the advertising content. In various embodiments, to accurately measure the influence of the advertising content data, the second user profile may be selected such that any attributes, characteristics, biases, confounding variables, and/or other factors that may affect the outcome of the measurement are reduced and/or eliminated. In certain cases, any variables potentially affecting the outcome of the measurement may be reduced by selecting a second user profile that is substantially similar (e.g., as close as possible) to the first user profile.

By way of example, a first user profile may include attribute data including demographic data (e.g., data indicating that the user is female, 30-40 years old, resides in San Francisco, Calif., has a household income of $100,000, etc.), behavioral data (e.g., the user visits a coffee shop three times per week), third party data (e.g., purchased a condo for $200,000 in 2006), psycho-graphic data (e.g., leads a healthy lifestyle, likely to vote for a particular political party, etc.), and other attribute data. Based on the attribute data, a second user profile that matches (e.g., is substantially similar to) the first user profile may be selected. The second user profile may include similar (e.g., matching) attribute data including demographic data (e.g., user is female, 30-40 years old, residing in San Francisco, Calif., household income of $95,000, etc.), behavioral data (e.g., visits the coffee shop four times per week), and/or other attribute data.

In one example matching approach, the attribute data from user profiles may be used in a regression approach (e.g., logistic regression, linear regression, etc.) to generate a model (e.g., generalized linear model (GLM), log it model, discrete choice model, etc.). For example, a model (e.g., generalized linear model (GLM)) may represent a correlation between a dependent variable of whether or not a user has been exposed to advertising content and a set/vector of covariates including attribute data included in the user profiles. The model (e.g., generalized linear model (GLM)) may be used to generate propensity scores for each of the multiple profiles. In some embodiments, a propensity score associated with the first user profile (e.g., associated with a user who has seen an ad) may be used to identify/select a matching (e.g., most closely matching) second user profile (e.g., associated with a user who has not seen the ad). A variety of matching approaches including nearest neighbor, kernel, local linear, caliper, and/or other matching techniques may be used to match the first and second user profiles based, for example, on propensity scores.

In one embodiment, first behavior information may be determined based at least in part on an association between the first user profile and a location associated with advertising content data. In various embodiments, behavior information may include information associated with a user's presence at one or more locations. In some embodiments, a first behavior information may include a number of instances, a number of instances over a period of time, and/or a frequency/rate that a user associated with the first user profile has been determined to be present at the location (e.g., visited the location). For example, a user may be determined to be present at a location based on location data (e.g., latitude/longitude and/or other location identifying information) received from a mobile device associated with the user. In certain cases, the location data may be received in connection with an advertisement request, a WiFi login page, marketing opportunity within a mobile application, entering a geo-fence, a deal and/or opportunity associated with a mobile device, etc. In various embodiments, location data received from a user device may be mapped to one or more defined locations. And based on a mapping of location data to a location associated with advertising content data, a user may be determined to be present at that location. When a user is determined to be present at a location, a user profile associated with that user may be updated to include information (e.g., behavioral information) associated with the user's presence at the location. For example, the user profile may be updated to include the location, a time (e.g., time/day stamp) of presence, duration of presence (e.g., five minutes), and/or other information related to the user's presence at the location. This information may be used to determine behavior information associated with the user profile and the location.

According to various embodiments, behavior information may include a number of times that and/or frequency with which a user associated with a user profile has been present at a location prior to and/or after being exposed to a digital advertisement. For example, a user associated with a first user profile may receive a digital advertisement including advertising content data associated with a location at certain time (e.g., a time (t0), a date, etc.). The time at which a user is exposed to advertising content data may include an advertising exposure time (e.g., time of exposure). In various embodiments, a user may have been exposed to advertising content data multiple times and the advertising exposure time may include the time of first exposure, time of last exposure, an average/median time over a period of multiple exposures, and/or any other time.

In some embodiments, behavior information associated with a first user profile may include a number of times a first user visited the location over a period of time (e.g., one week, three days, etc.) prior to exposure to advertising content data (e.g., viewing an ad). The period prior to exposure may include, for example, a look-back period. The look-back period may include any period of time (e.g., a predefined period, arbitrary period, etc.). A number, frequency, and/or rate at which a user visits a location during the look-back period may include a natural visit frequency/rate. A natural visit rate may represent a rate at which a user visits a location in the absence of exposure to advertising content (e.g., of the user's own volition, uninfluenced by advertising content, etc.).

In various embodiments, behavior information associated with the first user profile may include a number of times the user visited the location over a period of time after the time of exposure to the advertising content data (e.g., viewing the ad). The period of time after advertising exposure may include a look-forward period, and the look-forward period may be selected/determined in a manner similar to the look-back period. In certain cases, the look-forward period, however, may be selected to be substantially different than the look-back period. In another example, behavior information may include a frequency (e.g., one time per day, three times per week, etc.) at which the user visited the location during the look-forward period after exposure to the advertising content.

In some embodiments, behavior information may include a difference between a natural visit rate (e.g., a number of times and/or frequency at which a user was at the location during a period of time (e.g., a look-back period) prior to exposure to the advertising content data) and a number of times and/or frequency at which the user was at the location during a period of time after exposure (e.g., a look-forward period). The first behavior information may, for example, include value(s) quantifying an increase, decrease, and/or lack of change of the first user's behavior relative to the location (e.g., presence at the location) prior to and after seeing an advertisement. In various embodiments, an increase in presence at a location after viewing advertising content may indicate that the advertising content was successful in influencing the behavior of the user.

In various embodiments, behavior information may be determined based on location data from multiple mobile devices. For example, a user may be present at a location on a first day as determined by location information from a first device. After the first day, the user may replace the first device with a second device. Subsequently the user may be determined to be present at the location based on location data from the second device. In this case, location information received from both devices may be included in a user profile for the user, and behavior information may be determined based on location data from both devices that is included in the user profile.

In one embodiment, second behavior information may be determined based at least in part on an association between the second user profile and the location. In various embodiments, the second behavior information may include a number of instances, a number of instances over a period of time, and/or a frequency that a user associated with the second user profile (e.g., a control group profile) has been determined to be present at the location (e.g., visited the location).

In various embodiments, the second behavior information may include a change, if any, between the second user's visit frequency over a period (e.g., a look-back period) prior to a point in time as compared with the second user's visit frequency over a period (e.g., a look-forward period) after the point in time. The point in time (e.g., a reference time) may include, for example, the time at which the first user was exposed to the advertising content, a time relative to the time at which the first user was exposed to the advertising content, an arbitrary time, a time selected to ensure a proper comparison with the first behavior information, and/or another time.

In one embodiment, an advertising effectiveness value (e.g., a value representing advertising effectiveness, advertising effectiveness indicator) may be generated based at least in part on the first behavior information and the second behavior information. In some embodiments, an advertising effectiveness value may include number(s), value(s), percentage(s), metric(s) (e.g., a return on investment (ROI) metric, key performance indicator (KPI)), and/or any other data. The advertising effectiveness value may represent a change in number of visits (e.g., increase/lift in foot traffic) to a location as a result of exposure to the advertising content data.

In various embodiments, an advertising effectiveness value may be calculated/generated based on the first and second behavior information. In some embodiments, the advertising effectiveness value may be generated based on a comparison between a change in behavior from a time (e.g., a first time, a series of times, etc.) a first user sees an ad relative to their natural visit rate and a change in behavior of a second user who did not see the ad at the same time (e.g., an absolute same time, relative same time, etc.). Stated another way, the advertising effectiveness value may be generated based on a comparison of the first behavior information associated with a first user who saw an ad related to a location and second behavior information associated with a second user who did not see the ad. As discussed above, the first behavior information may include a change in a first user's visit behavior after exposure to advertising content relative to their natural visit rate. In other words, the first behavior information may be calculated based on a comparison (e.g., difference, change, etc.) of a first user's visit frequency to a location over a period of time (e.g., a look-back period) prior to exposure to advertising content related to the location and the user's visit frequency over a period after exposure (e.g., a look-forward period) to the advertising content. A second behavior information may include a change in behavior of a second user, who was not exposed to advertising content, as measured by a comparison of the second user's visit frequency to the location over a period (e.g., look-back period) prior to a certain time (e.g., the time when the first user saw the ad, a time relative to the time the first user saw the ad, an arbitrary time, etc.) and the second user's visit frequency over a period (e.g., look-forward period) after that time. The comparison of the first behavior information and second behavior information may be used to generate an incremental lift (e.g., advertising effectiveness value, which can be positive, negative, and/or zero) associated with the advertising content.

By way of example, first behavior information may indicate that a first user visited a coffee shop four times in the two weeks (e.g., a look-back period) prior to exposure to an ad for the coffee shop (e.g., an ad for a free coffee at the shop displayed to the first user on their mobile device). This visit rate over the look-back period (four times in two weeks (i.e., two times per week)) may include a natural visit rate for the first user. The first behavior information may also indicate that the first user visited the coffee shop four times in the week following exposure to the advertisement (e.g., a look-forward period). A second user profile may be matched to the first user profile using the matching techniques discussed herein. The second user may be a user with similar attributes to the first user. Second behavior information may indicate that the second user visited the coffee shop three times over the two weeks (e.g., a look-back period) prior to a point in time (e.g., the time the first user was exposed to the ad, a reference time, etc.) and two times in the week after that point in time. The advertising effectiveness value may be calculated based on the first behavior information and second behavior information. In one example, the advertising effectiveness value may include a comparison between a change in the first user's visit frequency prior to and after ad exposure time (e.g., four visits per week during the look-forward period versus two visits per week during the look-back period or a change/increase of two visits per week) and a change in the second user's visit frequency prior to and after the point in time (e.g., two times per week during the look-forward period and 1.5 times per week during the look-back period or a change of 0.5 visits per week).

In various embodiments, the process of generating advertising effectiveness values may be repeated for multiple pairs of users (e.g., associated with a location). And the multiple advertising effectiveness values may be aggregated (e.g., summed up, added together) to generate an aggregate advertising effectiveness value as discussed in detail below. An aggregate advertising effectiveness value including one or more advertising effectiveness values may include a location conversion index (LCI). In various embodiments, a group of users may be selected to determine an effectiveness/influence of advertising content (e.g., in driving users to a retail location). The group of users may, for example, be related to the location in some way (e.g., each user may have visited the location over a period of time, the users may have similar demographic attributes, etc.). The group of users may be divided into subgroups including an exposed subgroup (e.g., test subgroup) of users that have been exposed to the advertising content data and control subgroup including users not exposed to the advertising content data. Using the techniques discussed herein user profiles from the exposed subgroup may be paired to user profiles from the control subgroup and/or a general population subgroup. And advertising effectiveness values may be generated for each pairing of users, and the advertising effectiveness values may be aggregated (e.g., summed up) to generate an aggregate advertising effectiveness value. In various embodiments, the process of generating advertising effectiveness values may be performed iteratively across many different user profiles.

In some embodiments, the process of generating advertising effectiveness values may be repeated for multiple types of advertising content. For example, advertising effectiveness values may be generated for multiple versions of advertising content data.

In various embodiments, a user profile may include, for example, demographic data (e.g., household income, residence, value of home(s), occupation, work location, age, gender), behavioral data, data from third party data sources (e.g., property records, social network profile information, etc.), mobile device data (e.g., a list of applications on a device), psycho-graphic data, location visit frequency patterns, shopping cart spend data (e.g., including similar products and/or categories of products), and/or any other data associated with a user.

In some embodiments, behavioral attributes may be derived, for example, from a user's past locations (e.g., location pattern(s)), prior actions, and/or other data. For example, a user (e.g., associated with user profile) may have been determined to be at a location based on a location data received, for example, along with a mobile advertising request (e.g., from the user's mobile device). The location data may be mapped to a business, place of interest, zip+4 code, and/or other location. The mapped location data may be used to update a location pattern in the user's profile. The location patterns, behavior attributes, and/or other location-related information may be included in a location graph in, for example, the user's profile.

In some embodiments, demographic, behavioral, and/or other attributes associated with the business, place of interest, etc. to which a user's location has been mapped may be included in a user profile associated with that user. For example, a business (e.g., location) may be associated with demographic, behavioral, and/or other attributes. And as a result of a user's detected presence at the business, behavioral and/or other attributes associated with the business may be attributed to the user (e.g., added to a user profile associated with the user). In certain cases, attributes added to a user profile may be confirmed to be correct or incorrect based on other information (e.g., attributes associated with other locations the same user has visited, information from third party data sources, a user's device, etc.).

In some embodiments, an advertising effectiveness platform/service residing on one or more servers generates advertising effectiveness values (e.g., advertising effective index(es), location conversion index(es)/values, etc.) based on information derived from one or more user profiles. The advertising effectiveness service may query, mine and/or otherwise process user profile information stored in the user profile data store. For example, user profile information may be selected from the user profile data store and behavior information may be determined based on the selected user profile information. Advertising effectiveness values (e.g., generated based on the behavior information) may be stored in an advertising effectiveness data store. In various embodiments, an advertising provider may use the advertising effectiveness service to measure the effectiveness (e.g., influence, value, ROI, etc.) of an advertising campaign.

In one embodiment, an advertiser, advertisement provider, advertisement platform, and/or other entity may seek to determine an effectiveness of an advertising campaign associated with a retail location (e.g., an advertisement associated with a retail location). A first user may be selected based on a determination that the first user has been served advertising content associated with the campaign, the first user has visited the location prior to being served advertising content, and/or other criteria. In various embodiments, attribute data associated with a first user (e.g., included in a first user profile) may be used to select a second user. For example, location attribute data associated with the first user may indicate that the first user is a female, age 20-30, and employed at a technology firm. The location attribute data may also indicate that the first user visited the retail location (e.g., a fashion retailer) four times in the month prior to viewing an advertisement for the retail location. This natural visit frequency prior to being served the advertising content may include normal visits, unaided visits, and/or other types of visits to the retail location. Based on the first user's attribute data, a second user may be selected. In various embodiments, the second user may be selected using attribute-based matching, propensity score matching, and/or other matching approaches. The second user may, for example, include a user most similar (e.g., in demographic, behavioral, and/or other attributes; propensity score; and/or other metrics) to the first user. The second user may be selected based on a determination that the second user has not been exposed to the advertising content associated with the retail location and/or any advertising content associated with the retail location. In this example, a second user who is a female, age 20-30, employed at a law firm and visits the retail location three times per month may be selected. Whereas, another user who is a male, aged 40-50, employed as a doctor, and visits the retail location two times per quarter may not be selected as a similar user. The user may, however, be selected as a randomly-selected user as discussed below.

In various embodiments, first behavior information may be determined. In certain cases, the first behavior information may represent a comparison of a number of visits prior to and after the first user has been exposed to the advertising content (e.g., has viewed an ad, is presumed to have viewed an) associated with the retail location. According to some embodiments, second behavior information may be determined. In certain cases, the second behavior information may represent a number of times the second user visits the retail location prior to and after a certain point in time (e.g., the time the first user was exposed to the advertising content, another time, etc.). In various embodiments, an advertising effectiveness value may be generated based on the first behavior information associated with the first user and the second behavior information associated with second user. In various embodiments, the advertising effectiveness value may quantify/represent the influence of the advertising content data associated with the location.

According to some embodiments, an advertising effectiveness value may be generated based on a comparison of behavior information associated with the first user and behavior information associated with a randomly-selected user (e.g., a user from the general population). In various embodiments, a randomly-selected user may be selected based on a determination that the user is associated with the location (e.g., has visited the location over a period of time). It may be determined, for example, that the user has visited (360) the retail location; however, demographic data associated with user may not be similar to the demographic data associated with the first user. In various embodiments, an advertising effectiveness value may be generated based the behavior information associated with the first user and behavior information associated with the randomly-selected user using the approaches discussed herein. Generating an advertising effectiveness value based on a comparison of the behavior information associated with the first user and a randomly-selected user may provide additional insight into the effectiveness/influence of an advertisement.

In one embodiment, it may be determined that a user profile includes an indication of exposure to advertising content data and/or engagement to/with advertising content data. For example, an indication of exposure to advertising content data may include a record indicating that a digital advertisement including advertising content data associated with a location has been presented to a user. The indication may be generated, for example, when a digital advertisement is output to a user on a device (e.g., a mobile device, computer, smart television, wearable computer, etc.). An indication of engagement to/with advertising content data may include a record indicating that a user has engaged with advertising content by, for example, clicking on an ad, expanding an ad, engaging with an via voice input, and/or other records. In various embodiments, an indication of exposure/engagement may be associated with user profile and not a specific device. For example, a device (e.g., a home computer) on which the user was presented advertising content data and/or interacted with advertising content may be different than a device detected to be at a location of interest. In some embodiments, an indication of exposure/engagement may be generated when it is determined that a user has viewed and/or is likely to have viewed an advertisement presented in a non-digital medium (e.g., a print ad, mailed advertisement, etc.).

In one embodiment, the user profile may be selected based on the determination that the user associated with the profile has been exposed/engaged (e.g., is presumed to have viewed) to and/or engaged with the advertising content data. In various embodiments, a first user profile (e.g., test user profile) may be selected as a test user profile (e.g., for comparison with a control user profile as discussed herein) based on the determination that the first user profile includes an indication of exposure/engagement to the advertising content data.

In one embodiment, a continuity factor associated with a user profile may be determined. In various embodiments, continuity factors associated with user profiles may be used to select statistically significant user profiles. A continuity factor may indicate whether and/or to what extent a user was an active user (e.g., active in the system) prior to the time at which advertising content is served and/or after the advertising content has been served. A continuity factor, in some embodiments, may include a heart-beat indicator associated with the user. For example, if a user is determined to have been an active user on three separate days in the week prior to being served an advertisement for a location and three separate days after viewing the advertisement, the continuity factor for that user may be determined to be three. In various embodiments, the period of time prior to ad exposure and after ad exposure may be selected based on various factors associated with the advertising effectiveness calculation. The periods of time may, for example, be provided via user interface and/or other console from an advertiser.

In various embodiments, a continuity factor for a user profile may be generated based on location data from multiple mobile devices. For example, a user may be present at a location on a first day as determined by location information from a first device. After the first day the user may replace the first device with a second device. Subsequently the user may be determined to be present at the location based on location data from the second device. In this case, location information received from both devices may be included in a user profile for the user, and a continuity factor may be generated from the location data from both devices.

In one embodiment, it may be determined whether a continuity factor is above a threshold. In various embodiments, a threshold continuity factor may be set to, for example, one or any other value. A continuity factor greater than or equal to a threshold (e.g., one) may indicate that a user has been an active user before and after being served advertising content. This may indicate that the user profile is viable to be used in the propensity score calculation. In some embodiments, a continuity factor below a threshold (e.g., one) may indicate that the user was not present in the system prior to being served the advertisement. A user profile associated with a continuity factor below a threshold (e.g., one) may not be viable to be used in the propensity score calculation for purposes of evaluating the influence/effectiveness of advertising content data. In this case, the user may not be selected and the process may end.

In one embodiment, a user profile associated with a continuity factor above a threshold may be selected. In various embodiments, a user profile associated with a continuity factor value above a threshold may be selected as a test user profile (e.g., first user profile).

In one embodiment, propensity scores may be generated based on attribute data included in one or more user profiles. In some embodiments, a propensity score may represent a conditional probability of assignment to a particular treatment (e.g., exposure to the advertising content) given a set (e.g., vector) of observed covariates (e.g., attribute data including, for example, demographic attributes, behavioral attributes, psycho-graphic data, etc.). For example, a propensity score may represent a conditional probability of exposure to advertising content given a vector of attribute data (e.g., demographic data, behavioral data, psycho-graphic data, location visit frequency patterns, shopping cart spend data (e.g., including similar products and/or categories of products)).

In various embodiments, a propensity score associated with a user profile may be calculated by regressing the variable of whether or not the user has been exposed to advertising content against the attribute data included in the user profile. Using regression and/or other statistical approaches a model (e.g., generalized linear model (GLM), discreet choice model, etc.) may be generated representing a correlation between a dependent variable of whether or not a user has been exposed to advertising content and a set/vector of covariates including attribute data in the user profiles. In various embodiments, attribute data may be selected for inclusion in the set/vector of covariates to adjust for natural visit patterns, seasonal visit patterns, events, and/or other factors associated with the location of interest. The model (e.g., generalized linear model (GLM)) may be used to generate propensity scores for each of the multiple profiles. In some embodiments, the propensity score calculation process may account/compensate for natural visit patterns, seasonal visit patterns, events, and/or other factors associated with the location by virtue of the attribute data included in the propensity score calculation. For example, matching user profiles based on propensity score may reduce bias resulting natural visit patterns, seasonal visit patterns, events, and/or other factors.

In one embodiment, a first propensity score associated with the first user profile (e.g., a user profile in an exposed group) may be compared to one or more propensity scores each associated with a user profile in a control group (e.g., a group of user profiles for users not exposed to the ad content). In various embodiments, a first propensity score associated with the first user profile (e.g., a test group user profile) may be compared to one or more propensity scores to determine matching (e.g., closest/best matching) propensity scores.

In one embodiment, it may be determined whether a first propensity score matches one or more propensity scores. In some embodiments, a first propensity score may be compared to one or more propensity scores to determine a most-closely matching propensity score. In certain embodiments, nearest neighbor, kernel, local linear, caliper, and/or other matching techniques may be used to match the first propensity score to one or more propensity scores. In various embodiments, the first propensity score may be iteratively compared to multiple propensity scores to identify a most-closely matching propensity score. For example, a first propensity score (e.g., associated with a first user profile) may include a scalar value of 0.7, and this score may be compared to multiple propensity scores (e.g., 0.72, 0.65, 0.6, etc.) each associated with a user profile. Based on this example comparison, the propensity score of 0.72 may be selected as a most closely matching propensity score. In the event no propensity score is determined to match the first propensity score, the process may end.

In some embodiments, propensity scores may be matched based on a threshold and/or limit. For example, a first propensity score may match a second propensity score if the difference between the two propensity scores is within a threshold. For example, a first propensity score associated with a first user profile may include a scalar value of 0.35 and a second propensity score may include a scalar value of 0.3 and a threshold difference may be defined as 0.1. Because this difference between the first propensity score (e.g., 0.35) and second propensity score (e.g., 0.3) is less than the threshold (e.g., 0.1), the second propensity score may be determined to match (e.g., potentially match) the first propensity score.

In one embodiment, user profiles may be selected based on the matching propensity scores. In various embodiments, based on the propensity score matching process, the first user profile (e.g., including an indication of exposure to the advertising content) may be matched to a second user profile, and this pair of profiles may be selected. Once selected, an advertising effectiveness value may be calculated for the pair of user profiles.

In an embodiment of a process of calculating behavior information, a first timeline depicts a first user's behavioral patterns relative to a location (e.g., a retail location, restaurant, etc.) over a period of time. Each observation of the user (e.g., point) on the timeline may represent a point in time at which the first user was observed at the location. As depicted in the first timeline, the first user may, for example, have been served advertising content (e.g., associated with the location) at an ad exposure time (e.g., time of ad exposure, to, etc.). In some embodiments, a look-back period may include a period prior to the ad exposure time. A look-forward period may include a period after the ad exposure time. In some embodiments, the look-forward period and look-back period may include equal or different lengths/durations of time.

In some embodiments, first behavior information (e.g., associated with a user profile) may include a comparison of a first user's natural visit rate and post-advertising exposure visit rate (e.g., after exposure to the advertising content) to the location. A natural visit rate may include a number/frequency of user visits to the location over the look-back period. A post-exposure visit rate may include a number/frequency of visits to the location over the look-forward period after exposure to the advertising content. The first behavior information may include a difference (if any) between the first user's post-exposure visit rate and the natural visit rate.

In various embodiments, a second timeline is shown depicting a second user's behavioral patterns relative to a same location over a period of time. The second user in this case may not have been exposed to advertising content related to the location. In some embodiments, a look-back period for the second user may include a period prior to a point in time (e.g., a reference time). A look-forward period may include a period after the point in time. In various embodiments, the point in time (e.g., reference time) may be equivalent to the advertising exposure time (e.g., the same absolute time) at which the first user was exposed to the advertising content, another time determined based on the first and/or second user profile attributes, an arbitrary time, and/or any other time.

In some embodiments, the look-back period associated with the second user may be related to the look-back period associated with the first user. In one example, the two periods may span equivalent period(s) of time, though not necessarily the exact same period(s). For example, the first look-back period may include a first week (e.g., the last Wednesday in December to the first Wednesday in January, etc.), and the second look-back period may include (e.g., the first Saturday in February to the second Saturday in February). In another example, the first look-back period and second look-back period may span periods of time of varying duration. In various embodiments, similar relations may exist between the first look-forward period and second look-forward period.

In various embodiments, the look-back period, look-back period, look-forward period, look-forward period may be determined/selected based on input from a user of the advertising effectiveness platform, attributes associated with the first/second user profiles, and/or other parameters. In some embodiments, the look-back periods, and/or look-forward periods, may be selected to account/adjust for natural visit patterns, seasonal visit patterns, events (e.g., weather events, a sale at the location, etc.) associated with the location, and/or other factors that may influence/affect/skew the calculation of the advertising effectiveness value.

By way of example with reference to the first user timeline, a first user may be observed (e.g., via a mobile device) at a restaurant three times during the look-back period (e.g., as indicated by the three points on the timeline during the look-back period). The look-back period may include a one-week period prior to an ad exposure time of Jan. 1, 2014. The first user may have been shown advertising content related to the restaurant at the advertising exposure time (e.g., Jan. 1, 2014). And during the look-forward period including the two-week period after Jan. 1, 2014, the first user may be observed at the restaurant eight times. In this case, the first behavior information may include a difference between the first user's frequency of visits to the location during the look-back period—three times per week—and the first user's visit frequency during the look-forward period—four times per week. The first behavior information may include, for example, an increase of one visit per week, a 33.3% increase in visits per week, etc.

As depicted, for example, in the second user timeline, a second user may be observed at the restaurant (e.g., the same restaurant) four times during a second look-back period—the one-week period prior to Feb. 1, 2014. The second user may also be observed at the restaurant five times during a second look-forward period—the two weeks after Feb. 1, 2014. In this case, the second behavior information may include a difference between the second user's visit frequency to the location during the first look-back period—four visits per week—and the second user's visit frequency to the location during the second look-forward period —six visits over two weeks. The second behavior information may include, for example, a decrease of one visit of per week, a 25% decrease in visits per week, etc. In this case, the change in visit behavior after the reference time is negative (e.g., indicating a decrease). In certain cases, this negative value may be assumed to be the result from random behavioral patterns of the second user, and may be changed to zero indicating no change in behavior.

According to some embodiments, an advertising effectiveness value may be calculated based on the first behavior information and second behavior information. In this case the advertising effectiveness value may include a comparison between the first behavior information—an increase in one visit per week by the first user—and the second behavior information—a decrease of one visit per week by the second user. In this case, the advertising effectiveness value may include and incremental difference (e.g., incremental lift) of two visits per week. This value may indicate that exposure/interaction with the advertising content resulted in an increase visit frequency of two visits per week.

In one embodiment, two or more advertising effectiveness values may be generated. In various embodiments, a group of users including similar attributes may be selected to determine an effectiveness/influence of advertising content (e.g., in driving users to a retail location). For example, an advertiser associated with a quick service restaurant (QSR) chain may seek to quantify the value of an adverting campaign in driving foot traffic a QSR location. A group of user profiles identified as regular QSR patrons (e.g., known to visit the QSR location twice per week) may be selected. Within this group an exposed subgroup (e.g., exposed audience) of user profiles that include an indication of exposure to the advertising content may be identified. And a non-exposed subgroup of user profiles may be identified. Advertising effectiveness values may be generated using the techniques discussed herein. For example, user profiles from the exposed subgroup may be paired to similar user profiles from the non-exposed group, behavior information may be determined (e.g., numbers/frequencies of visits to the QSR location before and/or after advertisement exposure), and advertising effectiveness values may be generated based on the behavior information.

In one embodiment, aggregate effectiveness value(s) may be generated. In various embodiments, multiple advertising effectiveness values may be summed, aggregated, added together and/or otherwise combined to generate an aggregate advertising effectiveness value (e.g., a location conversion index). In various embodiments, an aggregate effectiveness value may include an advertising effectiveness value that has been updated based on other advertising effectiveness values. For example, two advertising effectiveness values may be merged/combined to generate a single advertising effectiveness value.

In various embodiments, advertising effectiveness values associated with any number of user profiles may be aggregated to generate the aggregate advertising effectiveness value. An aggregate advertising effectiveness value may represent an increase, decrease, and/or lack of change in a number of visits to retail location as a result of advertising content served to a defined group of users over a period of time. Continuing with the above example, the advertising effectiveness values generated based on the comparisons of the user profiles in the exposed subgroup and the users in the non-exposed subgroup of regular QSR patrons may be aggregated. For example, advertising effectiveness values may be generated for each user in the exposed subgroup and these values may be aggregated to generate an aggregate advertising effectiveness value across the group of regular QSR patrons. In one example, the aggregate advertising effectiveness value may, for example, indicate that the advertising campaign resulted in an increase of two visits per week per user who received the advertisement. In another example, the aggregate advertising effectiveness value may indicate a 25% increase in foot traffic to the QSR location over a defined period of time (e.g., one week before ad exposure compared to one week after ad exposure).

In various embodiments, advertising effectiveness values generated based on a comparison of user profiles exposed to advertising content and randomly-selected user profiles (e.g., not exposed to the ad content) may be included in an aggregate effectiveness value. For example, advertising effectiveness values may be generated based on comparisons of user profiles included in the exposed subgroup of male frequent QSR patrons to randomly-selected user profiles (e.g., not necessarily male frequent QSR patrons). These advertising effectiveness values may be added to an aggregate advertising effectiveness value, but may, for example, be given less weight in the aggregation.

In one embodiment, an aggregate advertising effectiveness value may be adjusted. In various embodiments, an aggregate advertising effectiveness value may be scaled, normalized, and/or otherwise adjusted. For example, advertising effectiveness value(s) may be scaled to a value within a range of values (e.g., 0 to 100), percentage(s), and/or other value(s).

In various embodiments, advertising effectiveness values may include adjustments for natural visit patterns, seasonal visit patterns, events, and/or other factors as a result of the matching processes (e.g., propensity score matching), look-back period determinations, look-forward period determinations, and/or other processes discussed herein. In some embodiments, however, an aggregate advertising effectiveness value (e.g., generated based on one or more advertising effectiveness values) may be adjusted (e.g., post-calculation) to account for natural visit patterns, seasonal visit patterns, events (e.g., current events, weather, etc.), and/or other factors affecting visit rates to a location. For example, an aggregate advertising effectiveness value reflecting ad campaign-driven visits to a retail location may be reduced to account for an increase in natural visits over the holiday season.

In one embodiment, a digital advertisement associated with a location may be generated. In various embodiments, a digital advertisement may include a coupon, a banner advertisement, a pop-up advertisement, embedded advertisement, and/or other promotional content associated with a location (e.g., aimed at driving foot traffic to the location). For example, a digital advertisement may include a coupon for a 20% discount on the purchase of a cup of coffee at a coffee shop.

In one embodiment, advertising effectiveness value(s) may be used to select users to receive the digital advertisement. In various embodiments, advertising effectiveness values may be used to select a type of user that would be most receptive to (e.g., most likely influenced by) the digital advertisement. Continuing with the example, an advertising effectiveness value may have been previously generated indicating that a coupon for a free muffin at the coffee shop resulted in an increased visit frequency of one visit per month among males, between 20-30 years old, with a median salary of $50,000 per year. Another advertising effectiveness value may have been generated indicating that a coupon for a 15% discount on purchase of coffee resulted in an increased visit frequency of two visits per week among males, between 40-50 years old, who regularly attend sporting events. Based on these advertising effectiveness values, user profiles associated with males, between 40-50 years, who are likely to attend sporting events may be selected to receive the digital advertisement.

In one embodiment, a digital advertisement may be provided to mobile device(s) associated with the selected user profiles. In various embodiments, providing digital advertisement to users in a group known to respond favorably to similar advertisement content may increase the return on investment of a mobile advertising campaign.

Further examples and details can be found in U.S. patent application Ser. No. 14/295,067, filed Jun. 3, 2014 and entitled “Measuring Advertising Effectiveness”, the entire disclosures of which application is hereby incorporated herein by references.

Measure of Mobile Visits Lifts

In one embodiment, location graphs are used to organize location data and attributes associated with mobile devices and predetermined geographical regions that may be visited by the mobile devices. A computing process is provided to propagate the attributes via the location graphs such that the attributes of the mobile devices and/or the predetermined geographical regions can be inferred according to a statistical model in a computationally efficient way. The attributes are used to quantify the characteristics of mobile devices and identify similar mobile devices. The impact of an event/information to the visitation patterns of the mobile devices is statistically measured as the differences between mobile devices that have been exposed to the event/information and similar mobile devices that have not been exposed to the event/information. Computing methods are provided to evaluate the impact statistically and provide a measurement of the impact.

Some details of location graphs can be found in U.S. Pat. App. Pub. No. 2014/0012806, entitled “Location Graph Based Derivation of Attributes”, the entire disclosure of which is hereby incorporated herein by reference. Further details can be found in the section below entitled “LOCATION GRAPH”.

Some details of the computation of the impact of an event/information to visitation patterns of mobile devices can be found in U.S. Pat. App. Pub. No. 2015/0348095, entitled “Measuring Advertising Effectiveness”, U.S. Pat. No. 9,374,671, entitled “Systems and Methods to Track Regions Visited by Mobile Devices and Detect Changes in Location Patterns”, and U.S. patent application Ser. No. 15/165,983, filed on May 26, 2016 and entitled “Systems and Methods to Track Regions Visited by Mobile Devices and Detect Changes in Location Patterns Based on Integration of Data from Different Sources”, the entire disclosures of which patent documents are hereby incorporated herein by reference.

In general, the users having similar or same profiles, including demographic profile data, behavioral profile data, psycho-graphic profile data, purchase profile data, and/or location profile data, etc., can be grouped as similar users. The similar users can be organized into two groups. One of the two groups is provided with a predetermined content (e.g., an advertisement, an announcement, a notice, a TV program, a direct mail advertisement, etc.); and the other of the two group is not provided with the content. The subsequent user behaviors, such as location patterns, web visitation, service subscription, retail store visitation, etc., can be compared between the groups to identify the influence of the content. In some embodiments, the changes in the user behaviors of the same users before and after the time of the presentation of the content are identified; and the changes are compared between the groups to identify the influence of the content.

Location Graph

It is challenging to identify, in a computer system, neighboring locations in a computational efficient way for a large set of locations.

At least one embodiment disclosed herein provides an efficient method for a computing system to identify neighboring relations among a set of locations on the surface of the Earth. The method uses a grid system, such as hierarchical grid systems illustrated in FIGS. 2-11, or other grid systems (e.g., Military Grid Reference System (MGRS)) to map locations to cells that contain the respective locations and to identify neighboring locations and/or candidates for neighboring locations based on cells that contain the respective locations and the neighboring cells.

FIGS. 23-25 illustrate a system to organize location data via a grid system according to one embodiment.

FIG. 23 illustrates a set (401) of locations (e.g., 411-415, 421-427). One embodiment disclosed herein provides a computational efficient way to identify the neighboring locations among the set (401) of locations (e.g., 411-415, 421-427) using a grid of cells as illustrated in FIG. 24.

FIG. 24 illustrates an operation of mapping locations (e.g., 411-415, 421-427) to the cells (e.g., 431-439) in a grid system. In FIG. 24, each location (e.g., 411-415, 421-427) is mapped to a corresponding cell (e.g., 431-439) that contains the location.

For example, in response to a determination that the locations (411-415) are in the cell (431), a set of data is stored in a computing system to facilitate the look up of the specific locations (411-415) are located within the cell (431). Thus, from the identifier of the cell (431), the set of locations that are located within the cell (431) can be looked up from the set of data stored in the computing system.

For example, the identifiers of the locations (411-415) can be stored in association with the identifier of the cell (431) in a lookup table such that the table can be queried using the identifier of the cell (131) to return the identifiers of the locations (411-415). For example, the computing device can be configured to store an array with cell identifiers as the indices of the array, and the lists of the identifiers of the locations (e.g., 411-415) contained within the respective cells as the values of the array for the corresponding indices. Other data storage techniques can also be used to facilitate the look up the locations that have been determined to be within the cell (431).

Preferably, the mapping of a location (e.g., 411) to a cell (e.g., 431) is performed via the direct manipulation of the coordinates of the location (e.g., 411) (e.g., based on resolution of the grid), without using a stored data table. For example, the mapping as illustrated in FIG. 13 can be used to convert the coordinates efficiently to the identifier of a cell at a desired grid resolution without using floating point number computations, when a hierarchical grid system as discussed in connection with FIGS. 8-13 is used. Since the coordinates of the location (411) can be efficiently converted to the identifier of the cell (431) that contains the location (411), it is not necessary to store data for the look up of the particular cell (431) that contains the location (411), when the grid reference system discussed above is used.

In one embodiment, for each respective location (e.g., 411) in the set of locations (e.g., 411-415, 421-427), the coordinates of the respective location (e.g., 411) are combined to generate the identifier of the cell (e.g., 431) that contains the respective location (e.g., 411). A data point is then added to the cell-location data to allow the subsequent look up of the respective location (e.g., 411) from the identifier of the cell (e.g., 431).

After the cell-location data is stored to facilitate the look up of locations contained within respective cells, by using the identifier of the cell as the index or query criterion, neighboring locations (or candidates of neighboring locations) of a location can be looked up from the identifier of the cell that contains the location and the identifiers of the neighboring cells.

For example, to identify the neighboring locations of any locations (411-415) in the cell (431), the computing system computes the identifiers of the neighboring/surrounding cells (e.g., 433-439) and then looks up, using the identifiers of the neighboring cells (e.g., 433-439) in the cell-location data, the locations (e.g., 421-427) that are contained within the neighboring cells. The collection of locations looked up for being in the cell (431) and its neighboring cells (433-439) identifies the neighboring locations (or candidates of neighboring locations).

Because the distances between a location (e.g., 411) within the cell (431) to any location (e.g., 413-415, 421-427) in the collection is less than a first threshold corresponding to a resolution of the grid (e.g., twice the length of the diagonal line of a grid cell), the collection of locations (e.g., 413-415, 421-427) are all within the first threshold from the location (e.g., 411).

Further, any location that is within a second threshold away from the given location (e.g., 411) are necessarily within the cell (e.g., 431) and its neighboring cells (e.g., 433-439), where the second threshold corresponds to a resolution of the grid (e.g., the length or height of the grid cell). The collection of locations (e.g., 413-415, 421-427) looked up from the identifiers of the cells (e.g., 431-439) includes all locations that are no more than the second threshold away from the given location (e.g., 411).

In some instances, a predetermined distance threshold can be used to select more precisely, from the collection of locations in the set of neighboring cells (431-433), the neighboring locations of the given location (411). For example, when the distance between a candidate location and the given location is no more than the distance threshold, the candidate location is selected as a neighboring location; otherwise, the candidate location is determined to be not a neighboring location.

The identification of the candidates based on the look up of locations from neighboring cells (431-433) reduces the number of candidates and thus the computation load for the comparison to the distance threshold. Preferably, the resolution of the grid is selected according to the predetermined distance threshold (e.g., in the same order as the threshold) to minimize the candidates that are not neighboring cells and/or minimize the neighboring cells that are to be looked up for candidates.

Preferable, the cell identifiers of the neighboring/surrounding cells (e.g., 433-435) of a given cell (431) can be computed directly from the identifier of the given cell (431) and/or the coordinates of a location (e.g., 411) inside the given cell (431) (e.g., as discussed in connection with FIG. 13). Thus, there is no need to store data for the look up of neighboring cells.

After the identification of the neighboring locations, the computer system stores graph data representing the neighboring relations among the locations (e.g., 411-415, 421-427) in the set (401) of locations. The graph data includes nodes representing the locations (e.g., 411-415, 421-427) and edges representing the neighboring relation between the locations and/or the distance between the neighboring locations, as illustrated in FIG. 25.

For example, in FIG. 25, the locations (411-415) are found to be located within a threshold distance from each other and hence connected via edges (441-445) in the graph of locations, where each of the locations is represented as a node in the graph. For example, in FIG. 25, the distance between locations (411 and 425) is less than the threshold and thus linked via an edge (447) in the location graph.

The graph data representing the neighboring relations among the locations can be stored using various techniques, such as look up table, linked lists, arrays, etc. The graph data allows the look up of the neighboring locations in the set (401) for any given location (e.g., 411).

FIG. 26 illustrates a location data processing system to establish a graph of locations according to one embodiment.

In the system illustrated in FIG. 26, a set of algorithms (e.g., 451, 453) are used to map locations to cells that contain the locations respectively, and map each cell to its neighboring cells. Through the mapping of the locations to cells using an algorithm (e.g., 451), cell-location data (e.g., 455) can be established and stored to facilitate the look up of locations contained within any cells. For each respective location, the cell contains the respective location and the neighboring cells are identified via the set of algorithms (e.g., 451, 453) and then used in the cell-location data (e.g., 455) to look up the locations that are contained within these cells. The looked up locations are identified as the neighboring locations (or candidates for the neighboring locations) of the respective location; and location graph data (457) can then be established and stored to facilitate the look up of neighboring locations of any location, in a way as illustrated in FIGS. 23-25.

In one embodiment, the algorithms (e.g., 451, 453) are based on a grid reference system (e.g., as illustrated in FIGS. 2-13). As illustrated in FIG. 13, the coordinates of a given location can be mapped to an identifier of a cell at a desired resolution, by manipulations of longitude digits and latitude digits to generate the column identifier and row identifier, which are further combined as a cell identifier. From the cell identifier, the coordinates of the vertices of the cell can be identified. Further, the coordinates of the vertices, as well as the identifiers, of the cells of the surrounding the cell can be identified based on the cell resolution. The description of FIG. 13 provides as further details.

Thus, the mapping relations from locations (e.g., 411, . . . , 415) to cells (e.g., 431) can be established based on the algorithm (451) applied to the coordinates of the locations (e.g., as illustrated in FIG. 13). Once a location (e.g., 411) is mapped to a cell (e.g., 431) that contains the location (e.g., 411), a portion of the cell-location data (e.g., 455) can be stored to allow the look up the location (e.g., 411) as part of the locations that are contained within the cell (e.g., 431). After all of the locations in a given set (e.g., 401) are mapped to the cells to store the relevant portions of the cell-location data (e.g., 455), the cell-location data (e.g., 455) can be used to look up all of the locations that are contained within any of the cells (e.g., 431).

Using the cell-location data (e.g., 455), the system then builds out the location graph data (e.g., 457) that connects any location (e.g., 411) to its neighboring locations (e.g., 413, . . . , 425). For example, for the location (411), the system identifies the cell (431) that contain the location (411) and the surrounding cells (e.g., 433-439), by applying the algorithms (e.g., 451, 453) to the coordinates of the location (411) and/or the identifier of the cell (431) in the grid reference system. The identifiers of the cell (431) and its neighboring/surrounding cells (e.g., 433-439) are used in the cell-location data (e.g., 455) to look up all of the locations contained within the cell (431) and its neighboring/surrounding cells (e.g., 433-439) for the identification of the neighboring locations (e.g., 413, . . . , 425) of the location (411).

Thus, in the system illustrated in FIG. 26, no data is required to be pre-stored for the mapping from the locations (e.g., 411-415) to the cells (e.g., 431) that contains the locations (e.g., 411-415); and no data is required to be pre-stored for the mapping from a cell (431) to its neighboring/surrounding cells (e.g., 433, . . . , 439).

For a given set (401) of locations (e.g., 411-415, 421-427) that are specified by their coordinates, the algorithm (e.g., 451) converts their coordinates to, in the grid reference system, the identifiers of the cells (e.g., 431) that contain the respective locations (e.g., 411-415). In response to the determination of the identifiers of the cells (e.g., 431) that contain the respective locations (e.g., 411-415), the cell-location data (e.g., 457) is stored to provide the mapping from the cells (e.g., 431) to the respective locations (411-415).

From the cell-location data that maps cells (e.g., 431) to respective locations (e.g., 411-415) that are contained within the respective cells (e.g., 431), the system identifies, for a given location (e.g., 411), the collection of neighboring locations (e.g., 413, . . . , 425) (or candidates of neighboring locations) that are located with a cell (431) that contains the given location (e.g., 411) and the neighboring/surrounding cells (433, . . . , 439). Optionally, the distances between the given location and the locations in the collection of candidates are computed and compared to a threshold to identify the neighboring locations (e.g., 413, . . . , 425).

Once the neighboring locations (e.g., 413, . . . , 425) of the given location (411) are identified, graph data (e.g., 457) is stored to map the location (411) to its neighboring locations (e.g., 413, . . . , 425). The process can be repeated for other locations (e.g., 415) to expend the graph data to include the mapping of any location to its neighboring locations.

FIG. 27 shows a method to generate a location graph according to one embodiment. For example, the method of FIG. 27 can be applied in the system of FIG. 1 using a grid reference system illustrated in FIGS. 2-13 to generate a location graph as illustrated in FIGS. 23-26.

In FIG. 27, a computing apparatus is configured to: receive (461) a set (401) of locations (e.g., 411-415, 421-427) each represented by its coordinates; compute, (463) for each location (e.g., 411) using its coordinates, the identifier of a cell (e.g., 431) that contains the location (e.g., 411) in a grid reference system (e.g., illustrated in FIGS. 8-13) and store cell-location data (e.g., 455) identifying, for each respective cell identifier of a cell that contains one or more of the locations in the set (401), the locations that are contained within the cell (e.g., 431) having the respective cell identifier; compute (465), for each respective location (e.g., 411), the cell identifier from its coordinates and identifiers of neighboring cells to look up, from the cell-location data (e.g., 455), locations identified by the computed cell identifiers; and compute (467) the distances between the respective location (e.g., 431) and the looked up locations (e.g., 413-415, 421-429) and store location graph data (e.g., 457) linking the respective location (411) with the neighboring locations (e.g., 413, . . . , 425) edges identifying the computed distances.

In one embodiment, all of the looked up locations (e.g., 413-415, 421-429) (other than the location (411) itself) are identified as neighboring locations. Alternatively, all of the looked up locations (e.g., 413-415, 421-429) (other than the location (411) itself) are identified as candidates of neighboring locations; and the distances between the location (411) and the candidates (e.g., 413-415, 421-429) are compared with a threshold to identify the neighboring locations (e.g., 413, . . . , 425) having distances to the location (411) that are no more than (or less than) the threshold.

The location graph data (457) can be used to propagate the attributes of the locations based on the proximity of the locations. For example, each of the locations in the set (401) may represent a business or a point of interest. The computing device may store a profile for each of the location in the set (401), where the profile identifies one or more known attributes (e.g., keywords) that persons visiting the location are likely to have. Due to the proximity of the locations, a person visiting a location is likely to visit a neighboring location; and thus a neighboring location is likely to have similar attributes. The location graph data (457) allows the system to propagate the profile attributes from locations to their neighboring locations based on the distances. For example, the likelihood of an attribute attachable to a person visiting the location can be computed based on the weighted average of the likelihood of the attribute in the profile of the location and the likelihoods of the attribute in the profiles of the neighboring locations identified in the location graph data (457), where the weights are based on the distances from the location to the neighboring locations. The longer the distance, the smaller the weight. The weighted average can be used to update the likelihood of the attribute in the profile of the location. The updating of the profile causes the propagation of the profile attributes from locations to neighboring locations.

The location graph can also be used to organize the locations tracked for a mobile device (e.g., tracked using a location determination device, such as a GPS receiver, of the mobile device).

In one embodiment, a computing device is configured (e.g., via software and/or hardware) to perform the operations to identify neighboring locations among a set of locations and/or create and store the location graph data.

For example, the computing device of one embodiment includes: at least one microprocessor; and memory storing instructions configured to instruct the at least one microprocessor to: store, in the computing device, coordinates of a plurality of locations (e.g., 411-415, 421-427) on a surface of the Earth. In a grid reference system, the surface of the Earth is covered by a plurality of cells (e.g., as illustrated in FIGS. 8-11. For each respective location in the plurality of locations, the instructions, when executed, cause the at least one microprocessor to: combine coordinates of the respective location into an identifier of a cell among the plurality of cells, where the cell contains the respective location on the surface of the Earth; and store, in the computing device, data associating the identifier of the cell and the respective location to facilitate a look up of the respective location using the identifier of the cell. For each respective location in the plurality of locations, the instructions, when executed, cause the at least one microprocessor to: identify a plurality of neighboring cells of the cell that contains the respective location on the surface of the Earth; look up, by the computing device, a subset of locations by using the identifier of the cell and the identifiers of the neighboring cells in stored cell-location data that associates identifiers of respective cells and locations contained within the respective cells; compute, by the computing device, distances between the respective location and locations in the subset; and store, in the computing device, graph data linking the respective location to locations in the subset with edges representing the distances, wherein the plurality of locations are represented as nodes in the graph data.

For example, the coordinates of the respective location are combined to provide the identifier of the cell according to a predetermined function of the coordinates of the respective location, without using additional data stored in the computing device, and without using a floating point number computation.

For example, the coordinates of the respective location are combined via: generating two integers from longitude and latitude coordinates of the respective location according to a precision level; and combining the two integers to provide the identifier of the cell.

For example, the coordinates of the respective location are combined via: selecting digits from the longitude and the latitude of the location in accordance with a cell resolution level; and combining the digits selected from the longitude and the latitude of the location into an integer representing the identifier of the cell.

For example, the plurality of neighboring cells are identified via computing identifiers of the neighboring cells from the identifier of the cell or the coordinates of the respective location; and the computing of the identifiers of the neighboring cells is based on predetermined functions of the identifier of the cell, without using additional data stored in the computing device and without a floating point number computation.

In one embodiment, a method implemented in the computing device, includes: storing, in the computing device, coordinates of a plurality of locations on a surface of the Earth, wherein the surface of the Earth is covered by a plurality of cells. For each respective location in the plurality of locations, the method further includes: combining, by the computing device, coordinates of the respective location into an identifier of a cell among the plurality of cells, wherein the cell contains the respective location on the surface of the Earth; and storing, in the computing device, data associating the identifier of the cell and the respective location to facilitate a look up of the respective location using the identifier of the cell. For each respective location in the plurality of locations, the method further includes: identifying, by the computing device, a plurality of neighboring cells of the cell that contains the respective location on the surface of the Earth; looking up, by the computing device, a subset of locations by using the identifier of the cell and the identifiers of the neighboring cells in stored cell-location data that associates identifiers of respective cells and locations contained within the respective cells; computing, by the computing device, distances between the respective location and locations in the subset; and storing, in the computing device, graph data linking the respective location to locations in the subset with edges representing the distances, wherein the plurality of locations are represented as nodes in the graph data.

In the method of one embodiment, the coordinates of the respective location are combined to provide the identifier of the cell without using additional data stored in the computing device.

In the method of one embodiment, the coordinates of the respective location are combined to provide the identifier of the cell according to a predetermined function of the coordinates of the respective location.

In the method of one embodiment, the combining of the coordinates of the respective location is performed by: generating two integers from longitude and latitude coordinates of the respective location according to a precision level; and combining the two integers to provide the identifier of the cell. For example, the two integers are combined to form the identifier of the cell without using a floating point number computation. In one embodiment, the cell is a rectangle area in a longitude latitude representation of the surface of the Earth.

In the method of one embodiment, the combining of the coordinates of the respective location is performed by: selecting digits from the longitude and the latitude of the location in accordance with a cell resolution level; and combining the digits selected from the longitude and the latitude of the location into an integer representing the identifier of the cell.

In the method of one embodiment, the identifying of the plurality of neighboring cells is performed by computing identifiers of the neighboring cells from the identifier of the cell. For example, the computing of the identifiers of the neighboring cells is based on the identifier of the cell without using additional data stored in the computing device. For example, the computing of the identifiers of the neighboring cells is based on predetermined functions of the identifier of the cell. For example, the computing of the identifiers of the neighboring cells is based on the identifier of the cell without a floating point number computation.

In the method of one embodiment, the identifying of the plurality of neighboring cells comprises computing identifiers of the neighboring cells from the coordinates of the respective location.

The method of one embodiment further includes: storing a set of keywords in association each of the plurality of locations; and propagating keywords associated with the plurality of locations via the edges in the graph data. For example, the propagating is performed based on weights computed based on distances represented by the edges in the graph data.

In one embodiment, a non-transitory computer storage medium storing instructions configured to instruct a computing device to perform any of the methods discussed above.

Keyword

Keywords can be used to describe the attributes of users of mobile devices (e.g., 109) and the attributes of users visiting predefined regions (e.g., 101, 103, . . . 105, 107). Examples of attributes include demographic attributes, behavioral attributes, propensity attributes, and/or psycho-graphic attributes.

For example, a set of keywords can be associated with the mobile device (109) as a profile of the user of the mobile device (109), which profile can also be considered as the profile of the mobile device (109). Similarly, a set of keywords can be associated with a predefined region (101) as a profile of users who visit the region (101), which profile can also be considered as the profile of the region (101).

Instead of simply indicating whether or not a profile has a particular keyword (e.g., a profile of a mobile device (109), or a profile of a region (101)), a profile can be configured to store a value for the particular keyword. The value or keyword value stored in the profile for a keyword indicates a level of confidence, probability, or likelihood that the profile has the attribute identified by the keyword. Such an arrangement allows the system to incrementally improve the accuracy of the profile using the methods discussed below.

For example, a value of zero for a keyword can be used in a profile to indicate that the likelihood of the profile having the keyword is zero. The keyword values can be scaled in a way such that all possible values are within a predetermined value range, such as the value range between zero and one, with a value of one for a keyword indicating that the probability of the profile having the keyword is one. Other value ranges can also be used based on suitable mappings or scales. A profile may store an ordered list of values corresponding to a predetermined list of keywords, or store keyword-value pairs for keywords that have values above a threshold (e.g., zero).

In one embodiment, the locations of the predefined regions and/or the mobile devices are organized into a location graph (e.g., using the technique discussed in connection with FIGS. 23-27). The location graph identifies, for each location (e.g., 411), its neighboring locations (e.g., 413, . . . , 425) and the distances from each location (e.g., 411) to its neighbors (e.g., 413, . . . , 425).

For example, a predefined region can be represented in a location graph by a location of its center. Alternatively, a predefined region can be represented in a location graph by the centers of cells of the grid reference system that contain portions of the predefined region. In one embodiment, each location (e.g., the coordinates of a location) is represented by an identifier of a cell at a suitable resolution of a grid reference system; and each predefined region is represented by the identifiers of a set of cells in the grid reference system, where the set of cells approximate the entire predefined region.

The locations of mobile devices (e.g., 109) (e.g., as identified by the coordinates of the mobile devices) can be mapped to the predefined regions (e.g., 101, 103, 105, 107) using the techniques such as those discussed in connection with FIGS. 12, 13, and/or 14. The location mapping technique allows the efficient counting of visits of each mobile device to the predetermined regions within a period of time (e.g., a day, a week). The visitation counts provide measurements of affinity of users of the mobile devices to respective regions.

After the neighbor relations among the predefined regions and distances to neighbor regions are computed for the location graph and the affinity of the users of the mobile devices to the respective regions (e.g., the form of visitation frequencies or visitation counts in a period of time), the inference of keyword values can be performed based on the neighboring relations, distances, and/or measured location affinity, as discussed below.

FIGS. 28 and 29 illustrate a method to propagate keyword values from neighboring locations based on distances to neighboring locations according to one embodiment.

FIG. 28 illustrates an example of a portion of a location graph that shows a location (501) and its neighbors (511, 512, 513, . . . , 519). The neighbor connections/links (e.g., 521, 522, 523, . . . , 529) indicate that the locations (511, 512, 513, . . . , 519) are the neighbors of the location (501), or neighbors of the first degree. Neighbors' neighbors are considered the neighbors of the second degree. A degree of separate between two locations in the location graph corresponds to the minimum number of neighbor connections/links in the location graph that connect the two locations. For example, if there are a minimum of n neighbor connections in the location graph to connect two locations, the two locations are separated by n degrees in the location graph.

In FIG. 28, each of the locations (501, 511, 512, 513, . . . , 519) has a keyword value (e.g., 502, 531, 532, 533, . . . , 539) for the same keyword (e.g., representative an attribute of users visiting the locations).

In FIGS. 28 and 29, the neighbor connections (521, 522, 523, . . . , 329) from the neighbors (511, 512, 513, . . . , 519) to the location (501) cause the neighbors (511, 512, 513, . . . , 519) to contribute keyword value updates for the location (501). The contributions are inversely proportional to the distances from the location (501) to the respective neighbors (511, 512, 513, . . . , 519). As illustrated in FIG. 29, the keyword value of the location (501) can be updated from an original value (502) to an updated value (503), which is a weighted average of the original value (502) and the values (531, 532, 533, . . . , 539) of the neighbors (511, 512, 513, . . . , 519) in the location graph.

For example, the weights (e.g., 541, 542, 543, . . . , 549) of the values (531, 532, 533, . . . , 539) of the neighbors (511, 512, 513, . . . , 519) can be evaluated according to a function of the respective distances of the neighbor connections (521, 522, 523, . . . , 529).

For example, a weight of one embodiment for a neighbor having a distance Di to the location (501) is proportional to 1/(1+Di/d) (and thus inversely proportional to Di), where d is a distance, such as the average distance from the neighbors (511, 512, 513, . . . , 519) to the location (501), the maximum (minimum, or medium) distance from the neighbors (511, 512, 513, . . . , 519) to the location (501), the threshold for determining whether a location is a neighbor, the maximum distance to a possible neighbor, the grid size used to generate the location graph, etc.

For example, a weight of another embodiment for a neighbor having a distance Di to the location (501) is proportional to 1/(b+Di/d), where d is a distance, similar to that discussed above; and b is a predetermined number. Thus, a closer neighbor of a location has a stronger influence/contribution in updating the keyword value of the location.

The sum of the weights of one embodiment equals to one in the weighted average (503) illustrated in FIG. 29. The unit of weight is split between the original value (502) and the contributions from the neighbors (511, 512, 513, . . . , 519) according to a predetermined ratio. For example, the sum of the weights (541, 542, 543, . . . , 549) for the values (531, 532, 533, . . . , 539) of the neighbors (511, 512, 153, . . . , 519) is equal to a predetermined number p, and the weight for the original value (502) of the location (501) is (1−p), where p is a number between 0 and 1, such as 0.5. The weights (541, 542, 543, . . . , 549) computed proportionally to the function of the distances (521, 522, 523, . . . , 529) can be scaled to have a sum that equals p, so that the sum of all weights in the weighted average (503) is one.

FIGS. 28 and 29 illustrate propagation of keyword values (531, 532, 533, . . . , 539) from neighbors (511, 512, 513, . . . , 519) of a location (501) according to the distances of the neighbor connections (521, 522, 523, . . . , 529) connecting the neighbors (511, 512, 513, . . . , 519) to the location (501). When the updating process is repeated x times, the contributions of neighbors having x degrees of neighbor separations can reach the location (501).

Further, the update (503) as illustrated in FIG. 29 can be optionally expanded to directly include the contribution of the neighbors of second degrees (e.g., neighbors of neighbors), where the sum of the weights for the neighbors of a high degree is scaled exponentially according to the degrees of separation. For example, the sum of the weights for the neighbors of x degrees can be p{circumflex over ( )}x (p to the power x). When the contributions of up to x degrees of neighbors are considered, the weight for the original value (502) of the location (501) is reduced to (1−p−p{circumflex over ( )}2− . . . −p{circumflex over ( )}x). Thus, the influence of the high order neighbors can be considered in the same iteration/update, with exponentially reduced influence according to the degrees of separations in the location graph.

Further, to reduce computation different neighbors of a same high degree can be given the same weight so that it is not necessary to compute the distances from the location (501) to its high order neighbors (e.g., neighbors of neighbors).

FIGS. 30 and 31 illustrate a method to propagate keyword values from neighboring locations based on visitation counts of locations according to one embodiment.

In FIGS. 30 and 31, the contributions of the values (531, 532, 533, . . . , 539) of the neighbors corresponding to a keyword to update the value (502) of the location (501) for the keyword are weighted based on the visitation counts (551, 552, 553, . . . , 559) of the neighbors (531, 532, 533, . . . , 539).

A visitation count for a region represents a measurement of intensity of mobile devices visiting the region. Multiple reports of a same mobile device to the region within a threshold period of time can be counted as a single visit.

In FIG. 31, the weights for the values (502, 531, 532, 533, . . . , 539) are directly proportional to the visitation counts (505, 551, 552, 553, . . . , 559) of the respective regions (501, 511, 512, 513, . . . , 519) and the sum of the weights equals to one.

The weighted average (507) illustrated in FIG. 31 can be further combined with the original value (502) according to a predetermined ratio. For example, the weighted average (507) is further provided with a weight of q, and the original value (502) with a weight of (1−q) for a weighted average as the updated value for the location (501), where q is a number between 0 and 1, such as 0.5.

The contributions of higher order neighbors can be considered in a way similar to that discussed above in connection with FIGS. 28 and 29 (e.g., through iterations or exponentially reduced weights according to degrees of separation).

FIGS. 32 and 33 illustrate a method to propagate keyword values from mobile devices to a location visited by the mobile devices according to one embodiment.

In FIGS. 32 and 33, the mobile devices (561, 562, 563, . . . , 569) are determined to have visited the location (501) with respective visitation counts (571, 572, 573, . . . , 579). The mobile devices (561, 562, 563, . . . , 569) have corresponding values (581, 582, 583, . . . , 589) for a keyword corresponding to the value (502) of the location (501). The values (581, 582, 583, . . . , 589) of the mobile devices (561, 562, 563, . . . , 569) are used to update the value (502) of the keyword for the location (501) based on their visitation counts (571, 572, 573, . . . , 579).

In FIG. 32, a count of a mobile device's visitation to the location (501) represents a measurement of intensity of the mobile device visiting the location (501). Multiple reports of the same mobile device to the region (501) within a threshold period of time can be counted as a single visit.

In FIG. 33, the weights for the values (581, 582, 583, . . . , 589) of the mobile devices (561, 562, 563, . . . , 569) are proportional to the visitation counts (571, 572, 573, . . . , 579) of the respective devices (561, 562, 563, . . . , 569) and the sum of the weights for the values (581, 582, 583, . . . , 589) of the mobile devices (561, 562, 563, . . . , 569) equals to a number r, which is a number between 0 and 1, such as 0.5. The weight for the original value (502) is (1−r), such that the sum of all of the weights in the update (509) equals to one. The update (509) can be used an updated value for the location (501).

The updates according to the methods of FIGS. 28-33 can be further combined. For example, the updates evaluated according to the method of FIGS. 32 and 33 for the neighbors (511, 512, 513, . . . , 519) and based on the values of mobile devices visiting the respective neighbors (511, 512, 513, . . . , 519) can be used as the values (531, 532, 533, . . . , 539) of the neighbors in computing the updates (503 or 507) according to the method of FIG. 29 or FIG. 31. For example, the updates (503, 507, and/or 509) computed according to the methods of FIG. 29, FIG. 31, and/or FIG. 33 (e.g., based on the values of the mobile devices having visited the location and/or the values of the neighbors of the location) can be combined via a weighted average according to a set of predetermined weights to update the value (502) of the location (501) for the keyword.

Thus, the methods discussed above allow the incremental improvement and enrichment of keyword values of the profiles of locations, based on inference from the attributes of neighbors and/or mobile device visitations to the locations.

FIGS. 34 and 35 illustrate a method to propagate keyword values from locations visited by a mobile device to the mobile device according to one embodiment.

In FIGS. 34 and 35, the value (601) of a keyword for a mobile device (109) is updated based on the values (621, 622, 623, . . . , 629) of the same keyword for regions/locations (611, 612, 613, . . . , 619), in accordance with the visitation counts (631, 632, 633, . . . , 639) of the mobile device (109) to the respective regions/locations (611, 612, 613, . . . , 619).

In FIG. 34, a count of visitations of the mobile device (109) to a region/location represents a measurement of intensity of the mobile device (109) visiting the region/location. Multiple reports of the same mobile device to the region (501) within a threshold period of time can be counted as a single visit.

In FIG. 35, the weights for the values (621, 622, 623, . . . , 629) of the regions/locations (611, 612, 613, . . . , 619) are proportional to the visitation counts (631, 632, 633, . . . , 639) of the device (109) to the respective regions/locations (611, 612, 613, . . . , 619); and the sum of the weights for the values (621, 622, 623, . . . , 629) of the regions/locations (611, 612, 613, . . . , 619) equals to a number s, which is a number between 0 and 1, such as 0.5. The weight for the original value (601) of the device (109) is (1−s), such that the sum of all of the weights in the update (603) equals to one.

Further, neighbors of the visited regions/locations (611, 612, 613, . . . , 619) (and/or their high order neighbors) can contribute to the values (621, 622, 623, . . . , 629) of the visited regions/locations (611, 612, 613, . . . , 619) (e.g., using the methods of FIGS. 29, 31 and/or 33) before the method of FIGS. 34 and 35 is applied to propagate the keyword values to the mobile device (109).

The method of FIGS. 34 and 35 allows the inference of the values of a mobile device (109) from the regions/locations visited by the mobile device (109); and the method of FIGS. 32 and 33 allows the inference of the values of a region/location from the mobile devices (e.g., 109) having visited the region/location. Further, the methods of FIGS. 28-31 allow the inference of values of a region/location from its neighbors in a location graph. In combination, the methods of FIGS. 28-35 allow the mutual improvements between the profiles of the regions/locations that indicate the attributes of users visiting the regions/locations and the profiles of users of mobile devices who visit various regions. The profiles can be further improved via information from other sources. When a profile is improved via another source (e.g., from a user input, or a third party data entry), the improvement is propagated via the visitation and location graph to locations visited, to neighbor locations, and to mobile devices visited the locations having the improved profiles).

FIG. 36 shows a method to update keyword values according to one embodiment.

In FIG. 36, a computing device (e.g., the server (187) illustrated in FIG. 1) is configured to: store (641) (e.g., in the database (181) illustrated in FIG. 1) a first set of location data identifying a plurality of predefined locations each having a predefined boundary; receive (643), from a plurality of mobile devices, a second set of location data identifying locations of the mobile devices during a period of time; store (645) (e.g., in the database (181) illustrated in FIG. 1), for each keyword identifying a user attribute, keyword values of the plurality predefined locations respectively and keyword values of the plurality of mobile devices respectively; generate (647) a location graph that has neighbor links connecting neighboring ones of the predefined locations and distances between locations directly connected by the neighbor links; identify (649), using a grid reference system, corresponding locations in the predefined locations that contain respectively the locations of the mobile device; determine (651), based on a result of identifying the corresponding locations, visitation counts of the mobile devices to the predefined locations during the period of time; and update (653), for each keyword, the keyword values of the predefined locations and the mobile devices based on the visitation counts and the location graph.

The location graph can be generated using the technique discussed in connection with FIGS. 23-27.

The identification of the respective predefined locations that contain the respective locations of the mobile devices can be performed using the technique discussed in connection with FIGS. 12-14.

The update (653) can be performed in a way as illustrated in FIGS. 28 and 29, FIGS. 30 and 31, FIGS. 32 and 33, FIGS. 34 and 35, or combinations.

For example, the plurality of predefined locations includes a first predefined location (501) and second predefined locations (511, 512, 513, . . . , 519) that are connected to the first predefined location directly via neighbor links (521, 522, 523, . . . , 529) in the location graph illustrated in FIG. 28 or 30.

To update a value (502) of the first predefined location (501) for the keyword based on values (531, 532, 533, . . . , 539) of the second predefined locations (511, 512, 513, . . . , 519) for the keyword, the computing device computes a weighted average of the values (531, 532, 533, . . . , 539) of the second predefined locations (511, 512, 513, . . . , 519), as illustrated in FIG. 29 or 31.

According to FIG. 29, the weights (541, 542, 543, . . . , 549) are computed computing, for the weighted average of the values (531, 532, 533, . . . , 539) of the second predefined locations (511, 512, 513, . . . , 519) based on distances (D1, D2, D3, . . . , Dm) from the first predefined location (501) to the second predefined locations (511, 512, 513, . . . , 519) respectively. The weights are inversely proportional to the distances (D1, D2, D3, . . . , Dm) from the first predefined location (501) to the second predefined locations (511, 512, 513, . . . , 519).

According to FIG. 31, the weights for the values (531, 532, 533, . . . , 539) of the second predefined locations (511, 512, 513, . . . , 519) are computed based on visitation counts (551, 552, 553, . . . , 559) of the mobile devices to the second predefined locations (511, 512, 513, . . . , 519) respectively. The weights are directly proportional to the visitation counts (551, 552, 553, . . . , 559) of the mobile devices to the second predefined locations (511, 512, 513, . . . , 519). Further, the weighted average can include a contribution from the value (502) of first predefined location (501) weighted according to a visitation count (505) of the mobile devices to the first predefined location (501).

The updates computed according to FIG. 29 and FIG. 31 can be combined via weighted average according to a set of predetermined weights.

Further, the keyword values of mobile devices that have visited the predefined locations can contribute to the update (6530, as illustrated in FIGS. 32 and 33. For example, the plurality of mobile devices include a subset of first mobile devices (561, 562, 563, . . . , 569) that have visited the first predefined location (501). To update the value (502) of the first predefined location (501) using the values (581, 582, 583, . . . , 589) of the first mobile devices (561, 562, 563, . . . , 569) for the keyword, the computing device computes a weighted average of the values (581, 582, 583, . . . , 589) of the first mobile devices (561, 562, 563, . . . , 569) as illustrated in FIG. 33. The weights for the values (581, 582, 583, . . . , 589) of the first mobile devices (561, 562, 563, . . . , 569) are based on visitation counts (571, 572, 573, . . . , 579) of the first mobile devices (561, 562, 563, . . . , 569) respectively to the first predefined location (501). For example, the weights for the values (581, 582, 583, . . . , 589) of the first mobile devices (561, 562, 563, . . . , 569) are computed to be directly proportional to the visitation counts (571, 572, 573, . . . , 579) of the first mobile devices (561, 562, 563, . . . , 569) to the first predefined location (501).

The update computed according to FIG. 33 can be combined with the updates from FIG. 29 and/or FIG. 31 through a weighted average according to a set of predetermined weights.

Further, the updates according to FIG. 33 can be performed for the second predefined locations (511, 512, 513, . . . , 519) before the update according to FIG. 29 and/or FIG. 31 is computed.

Furthermore, the keyword values (621, 622, 623, . . . , 629) of the predefined locations (611, 612, 613, . . . , 619) that have been visited by the mobile device (109) can be used to update the keyword value (601) of the mobile device (109) via a weighted average, as illustrated in FIGS. 34 and 35, where the weights for the predefined locations (611, 612, 613, . . . , 619) are directly proportional to the visitation counts (631, 632, 633, . . . , 639) of the mobile device (109) to the predefined locations (611, 612, 613, . . . , 619).

In one embodiment, each keyword value (e.g., 502, 531, 532, . . . , or 539) of a predefined region (e.g., 501, 511, 512, . . . , or 519) represents the likelihood (or probability or confidence level) of users who visit the respective region (e.g., 501, 511, 512, . . . , or 519) having the user attribute represented by the keyword. Each keyword value (e.g., 581, 582, . . . , or 589) of a mobile device (e.g., 561, 562, . . . , or 569) represents the likelihood (or probability or confidence level) of the user of the respective mobile device (e.g., 561, 562, . . . , or 569) having the user attribute represented by the keyword. Through an iterative process of propagating the keyword values from mobile devices and/or neighbor locations, and from regions to mobile devices that have visited the regions, the accuracy of the likelihood (or probability or confidence level) increases.

FIGS. 28-36 illustrate the propagation of the values for one keyword as an example. The values of a plurality of keywords can be propagated in a similar way.

Time Series

In one embodiment, time-stamped location information of mobile devices is processed to generate a statistical measurement of a change in location behaviors of mobile devices after a time period of influence related to a predetermined region.

For example, the server (187) is configured in one embodiment to receive location inputs regarding a population of mobile devices (e.g., 109 illustrated in FIG. 1). Each location input identifies a mobile device (e.g., 109), a location (111) of the mobile device (e.g., 109), and the date and time of the location (111). In general, the location (111) may or may not be inside the predetermined region (e.g., 101 or 103). The location input may be communicated from the mobile devices (e.g., 109) to the server (187) directly over the network (189). Alternatively, the location input may be communicated from the mobile devices (e.g., 109) to the server (187) indirectly via the database of a separate computing system (e.g., for enhanced privacy protection). An example of such a communication method can be found in U.S. patent application Ser. No. 15/165,983, filed May 26, 2016 and entitled “Systems and Methods to Track Regions Visited by Mobile Devices and Detect Changes in Location Patterns based on Integration of Data from Different Sources”, the entire disclosure of which is hereby incorporated herein by reference.

During a time period of influence, one subset of the population of the mobile devices (e.g., 109) is subjected to influences related to the predetermined region (e.g., 101), while another subset of the population of the mobile devices (e.g., 109) is not subjected to the influences related to the predetermined region (e.g., 101). The server (187) converts the location inputs into a statistical measurement that quantifies the effect of the influences on the location behaviors of the mobile devices (e.g., 109).

To quantify the influences, the server (187) converts the location inputs of a respective mobile device (e.g., 109) into measurements of visitations of the mobile device (e.g., 109) to the predetermined region (101) for a plurality of time intervals. Each of the time series of visitation measurements represents the intensity of visits to the predetermined region (101) within a time interval of a predetermined length (e.g., a hour, a day, a week, or a month). For example, the server (187) uses the method of FIG. 12 or 14 to determine whether a respective location of the mobile device (e.g., 109) is with the predetermined region (101), and combines a sequence of time instances of the mobile device (e.g., 109) being located within the predetermined region (101) during the time interval into a visitation measurement for the time interval (e.g., by combining time instances that are separated by a time period less than a threshold as one visit, and counting visits of the mobile device (109) in the time interval).

The server (187) applies a weight function to the time series of visitation measurements to compute a change of visitation intensity before and after a time instance. The change of visitation intensity before and after a time instance is computed from a weighted average of visitation measurements for one or more time intervals before the time instance and the visitation measurements for the same number of time intervals after the time instances. The weights computed according to the weight function for the time intervals of the visitation measurements decrease (e.g., exponentially) as the time gap between the time instance and the time intervals of the visitation measurements increases. The weights for the time intervals after the time instance has one sign (e.g., positive) and the weights for the time intervals before the time instance has the opposite sign (e.g., negative), such that the visitation measurements before and after the time instance subtracts from each other to indicate the change before and after the time instance.

The server (187) sums the visitation changes of the mobile device (109) computed for a plurality of time instances during the time period of influence as a response of the mobile device (109) for the time period. The server (187) computes an average of the responses of exposed mobile devices, which are exposed to the influences in the time period of influence, and an average of the responses of the non-expose mobile devices, which are not exposed to the influences in the time period of influence. The difference between the average responses of the exposed mobile devices and non-exposed mobile devices provides a measurement of a change in location behaviors that is the result of the influences applied during the time period of influence.

Preferably, the set of non-exposed mobile devices which are not exposed to the influences in the time period of influence are selected to be similar to the set of exposed mobile devices that are exposed to the influences in the time period of influence. The similarity can be identified based on the keyword values of the mobile devices before the time period of influence (e.g., the keyword values determined using the method of FIG. 36 and/or FIGS. 28-35).

For example, when a measurement of the differences between the keyword values of two mobile devices are within a threshold, the two mobile devices can be considered to be similar to each other. For example, a best matching mobile device can be found for a matched mobile device by searching for a mobile device that has the smallest measurement of the differences between the keyword values of the matching and matched mobile devices; and the best matching mobile device can be considered the mobile device that is most similar to the matched mobile device.

In one embodiment, for each exposed mobile device in the exposed set of mobile devices, the server (187) identifies, from non-exposed mobile devices, a non-exposed mobile device that is similar to (matching with, best matching with, or most similar to) the respective mobile device in the exposed set. The difference between the responses of the pair of exposed and non-exposed mobile devices is an indication of the behavior change resulting from the time period of influence; and an average of such differences for pairs of exposed and non-exposed mobile devices provides a statistical measurement of the behavior change, which is equal to the difference between the average response of the exposed mobile devices and the average response of the non-exposed mobile devices.

In general, the non-exposed mobile devices selected in the computation of the statistical measurement of a change in location behaviors do not have to match on a one-to-one basis with the exposed mobile devices used in the computation. For example, there may be more or fewer of the non-exposed mobile devices used in the computation of the average response of non-exposed mobile devices than the exposed mobile devices used in the computation of the average response of exposed mobile devices; and the difference between the average response of exposed mobile devices and the average response of non-exposed mobile devices provides a statistical measurement of the change in location behaviors. In some embodiments, the non-exposed mobile devices are selected to match the distribution of the exposed mobile devices in a space of keyword values such that a statistical distribution of the non-exposed mobile devices in the space of keyword values is the same as, or approximately the same as, the statistical distribution of the exposed mobile devices in the space of keyword values.

FIG. 37 illustrates a method to construct a visitation measurement for a time interval according to one embodiment. In FIG. 37, a time interval between time instances Ta and Tb are divided into a plurality of segments (e.g., t1, t2, . . . , tn) (e.g., based on a predetermined time period threshold). When a mobile device (109) has at least one location being within a predetermined region (101) during a time segment (e.g., t1, t3, . . . , or tn), the mobile device (109) is considered to have one visit (e.g., 671, 673, . . . , or 675) to the predetermined region (101) corresponding to the time segment (e.g., t1, t3, . . . , or tn). The total count of visits (671, 573, . . . , 675) for the time interval between time instances Ta and Tb provides a visitation measurement Vb (677) for the time interval. The visitation measurement Vb (677) corresponds to the total count of time segments (t1, t3, . . . , or tn) during each of which time segments the mobile device (109) has at least one location in the predetermined region (101). The visitation measurements can be scaled from the visitation counts so that the visitation measurements is within a predetermined data range (e.g., between 0 and 1).

Alternatively, for example, the server (186) generates, based on the location inputs of the mobile device (109), a sequence of the time instances of locations of the mobile device (109) that are located within the predetermined region (101) and fall between the time instances Ta and Tb. When two adjacent time instances have a time gap shorter than the predetermined time period threshold, one of the time instances is removed. The process is repeated until no two adjacent time instances have a time gap shorter than the predetermined time period threshold. The visitation measurement Vb is proportional to a count of remaining time instances in the sequence.

Alternatively, for example, the server (186) generates, based on the location inputs of the mobile device (109), a sequence of the time instances of locations of the mobile device (109) by adding time instances one at a time. When a new time instance of the mobile device (109) in the predetermined region (101) is within the predetermined time period threshold from an existing time instance in the sequence, the new time instance is not added to the sequence; otherwise, the new time instance is added to the sequence. The process is repeated until all location data of the mobile device (109) have been processed. The visitation measurement Vb is proportional to a count of time instances in the sequence that are between the time instances Ta and Tb.

Alternatively, for example, the server (186) generates, based on the location inputs of the mobile device (109), a sequence of the time instances of locations of the mobile device (109) that are located within the predetermined region (101) and fall in the time interval from Ta to Tb. The server (186) performs a statistical analysis of the sequence of the time instances to identify a plurality of clusters of time instances. The visitation measurement Vb is proportional to a count of clusters of time instances.

After using one of the methods to convert the location data of the mobile device (109) to a time series of visitation measurements for a plurality of time intervals of a predetermined length, the server (186) computes a response for a time instance based on one or more visitation measurements before the time instance and corresponding one or more visitation measurements after the time instance, as illustrated in FIG. 38.

FIG. 38 illustrates a method to construct a location response at a time instance from visitation measurements of a plurality of time interval according to one embodiment.

In FIG. 38, a weight function (683 and 685) is centered at a time instance (681) between the time intervals. The weight function (683 and 685) is anti-symmetric with respect to the time instance (681). For visitation measurements (e.g., 676 and 677) of time intervals that have equal elapsed times before and after the time instance (681), the weight function (683 and 685) provides weights (682 and 684) of the same magnitude but opposite signs. Thus, the weighted average to the respective visitation measurements (e.g., 676 and 677) represented a weighted difference between the respective visitation measurements (e.g., 676 and 677) that are have equal elapsed times before and after the time instance (681). The magnitude of the weight function (683 and 685) decreases (e.g., exponentially) as the elapsed time from the time instance (681) increases. Thus, the weighted average of visitation measurements computed according to the weight function correspond to the sum of the weighted differences in pairs (e.g., 676 and 677) of visitation measurements having equal elapsed times from the time instance (681), where the weights decrease as the elapsed times increase. The weighted average represents a location response (687) of the location behavior relative to the time instance (681).

The computation of the location response (687) of the time instance (681) can be repeated for other time instances between adjacent time intervals of the time series of visitation measurements, as illustrated in FIG. 39.

FIG. 39 illustrates a method to construct a location response for a time period according to one embodiment.

In FIG. 39, the mobile device (109) has a sequence of location responses (687, 688, . . . , 689) at respective time instances within the time period (691) of influence. Each of the location responses (687, 688, . . . , 689) is between two adjacent time intervals of the time series of visitation measurements. The server (187) computes a sum of the sequence of location responses (687, 688, . . . , 689) during the time period (691) of influence as the response (693) of the mobile device (109).

After computing the responses (e.g., 693) of the exposed mobile devices (109) and the corresponding responses of the non-exposed mobile devices, the server (187) generates a statistical measurement of a change in location behaviors from the difference between the average of the responses (e.g., 693) of the exposed mobile devices (109) and the average of the corresponding responses of the non-exposed mobile devices. In some instances, when there is a one-to-one mapping between the exposed mobile devices (e.g., 109) and the non-exposed mobile devices that are similar to the respective exposed mobile devices (e.g., 109) based on their keyword values, the server (187) generates the statistical measurement of the change in location behaviors from the average of the differences between the responses (e.g., 693) of the exposed mobile devices and the respective responses of the non-exposed mobile devices.

FIG. 40 shows a method to determine a measurement of a location change according to one embodiment. For example, the method of FIG. 40 can be implemented in a system illustrated in FIG. 1.

In FIG. 40, a server (e.g., 187) is configured to convert (701) a set of tracked locations of a mobile device (109) into a time series of measurements of visitations to a predetermined region (101). Each of the visitation measurements (e.g., 677) is computed for a time interval (e.g., Ta to Tb) using a method as illustrated in FIG. 37, or another method discussed above).

The server (e.g., 187) is further configured to apply (703) a weight function (e.g., 685 and 683) to a number of visitation measurements (e.g., 677) before a time instance (681) and the same number of visitation measurements (e.g., 676) after the time instance (681) to compute a response (687) of the mobile device (101) for the time instance (681), in a way as illustrated in FIG. 38.

The computation of a response at a time instance (681) is repeated for a plurality of time instances corresponding to a set of time intervals within a predetermined time period (691) of influence as illustrated in FIG. 39.

The server (e.g., 187) is further configured to sum (705) the responses (687, 688, . . . , 689) of the mobile device (109) for the plurality of time instances within the predetermined time period (691) as a measurement of a change response (693) in location behavior of the mobile device (109) resulting from one or more influencing events in the predetermined time period (691).

FIG. 41 shows a method to measure a statistical change in location behavior of mobile devices according to one embodiment. For example, the method of FIG. 41 can be implemented in a computing device, such as a server (187) illustrated in FIG. 1.

In FIG. 41, the computing device is configured to: receive (711) locations of mobile devices (e.g., 109) at different time instances; receive (713) data identifying a time period (691) of influence related to a predetermined region (e.g., 101); identify (715) exposed mobile devices that are exposed to the influences during the time period; identify (717) non-exposed mobile devices that are not exposed to the influences during the time period; for (719) each respective mobile device (e.g., 109) in the exposed mobile devices and the non-exposed devices, determine a set of time instances at which the respective mobile device (e.g., 109) has locations in the predetermined region (e.g., 101) at the respective time instances in the set; convert, (721) using an anti-symmetric weight function having an exponential distribution (683 and 685), the time instances into a response (693) in location behavior change of the respective mobile device (e.g., 109) for the time period (691) of influence; and compute (723) a statistical measurement of changes resulting from the time period (691) of influence from a difference between the average response of the exposed mobile devices and the average response of the non-exposed mobile devices.

For example, the system illustrated in FIG. 1 has a plurality of mobile devices (e.g., 109) and a computing device (e.g., the server (187)) in communication with the mobile devices (e.g. 109) to receive locations of the mobile devices at different time instances.

The computing device (e.g., 187) receives information identifying a time period (691) during which at least one message related to a predetermined region (e.g., 101) is transmitted to a portion of the mobile devices at least once to influence the location behaviors of the exposed mobile devices that have received the message.

The computing device (e.g., 187) identifies a first subset of the mobile devices, where each mobile device in the first subset is an exposed mobile device that has received the message at least once during the time period (691).

The computing device (e.g., 187) also identifies a second subset of the mobile devices, where each mobile device in the second subset is an non-exposed mobile device that has not received the message during the time period. The second subset can be identified based on similarity to the first subset. For example, for each exposed mobile device in the first subset, a non-exposed mobile device having keyword values closest to those of the exposed mobile device at a time before the time period can be selected and added to the second subset, such that there is a one-to-one mapping between exposed mobile devices in the first subset and the closest matching non-exposed mobile devices the second subset. Alternatively, a non-exposed mobile device is selected to have a distribution in a keyword value space that best matches with the distribution of the exposed mobile device, without requiring a one-to-one mapping between the first subset and the second subset.

For each respective mobile device (e.g., 109) in the first subset and the second subset, the computing device (e.g., 187) determines whether or not each location, in a plurality of locations of the respective mobile device (109) at a plurality of time instances respectively, is located within the predetermined region. The determination can be performed efficiently using the method of FIG. 14 or FIG. 16.

The computing device (e.g., 187) then converts the plurality of locations of the respective mobile device (e.g., 109) into a time series of measurements (e.g., 677) of visitation of the respective mobile device (109) to the predetermined region (101) for a plurality of time intervals (e.g., Ta to Tb) corresponding to the time series.

The computing device (e.g., 187) computes a plurality of responses (e.g., 687) at a plurality of time instances (e.g., 681) during the time period (691) by applying a weight function (683 and 685) to the time series of the measurements (e.g., 677, 676) of visitation to the predetermined region (101).

After the computing device (e.g., 187) sums the plurality of responses (e.g., 687, 688, . . . , 689) at the plurality of time instances as a response (693) of the respective mobile device (109), the computing device (e.g., 187) determines a statistical measurement of a location behavior change based on a difference between: an average of responses of mobile devices in the first subset; and an average of responses of mobile devices in the second subset.

The computing device (e.g., 187) of one embodiment converts the plurality of locations of the respective mobile device (109) into the time series of measurements of visitation to the predetermined region by counting multiple locations of the respective device located within the predetermined region and within a predetermined time period as one visit to the predetermined region.

For example, each of the measurements of visitation to the predetermined region in the time series can be based on a count of visits to the predetermined region within a respective time interval of the plurality of time intervals. The measurements of visitation to the predetermined region are the counts of visits normalized to have values that are within a predetermined range (e.g., between 0 and 1).

For example, each measurement, in the time series of measurements of visitation to the predetermined region, computed for a respective time interval of the plurality of time intervals, is proportional to a count of visits of the respective mobile device to the predetermined region during the respective time interval.

In one embodiment, each respective time interval of the plurality of time intervals has a plurality of time segments of a predetermined length; and each measurement, in the time series of measurements of visitation to the predetermined region, computed for the respective time interval, is proportional to a count of time segments, in the respective time interval, during which time segments the respective mobile device has at least one location that is located within the predetermined region.

Alternatively, each measurement, in the time series of measurements of visitation to the predetermined region, computed for a respective time interval of the plurality of time intervals, is proportional to a count of statistical clusters of locations of the respective mobile device in the respective time interview, wherein the clusters are at least a predetermined time duration apart in time. For example, the computing device of one embodiment performs a statistical cluster analysis of locations of the respective mobile device to identify the statistical clusters.

To apply the weight function (683, 685) to the time series of the measurements (677, 676) of visitation to the predetermined region (101), the computing device computes, for each respective time instance (e.g., 681) in the plurality of time instances, a weighted average of a subset of the time series of the measurements according to the weight function. The subset of the time series of the measurements includes a number of first measurements (e.g., 677) for time intervals that are within a predetermined time length (e.g., the time span of the portion (683) of the weight function) before the respective time instance (681) and the same number of second measurements for time intervals that are within the predetermined time length (e.g., the time span of the portion (685) of the weight function) after the respective time instance (681). The weights for the first measurements (e.g., 677) and weights for the second measurements (e.g., 676) are computed from the weight function (683 and 685) based on elapsed times between the respective time instance (681) and time intervals of respective ones (e.g., 677 and 676) of the first measurements (e.g., 677) and the second measurements (e.g., 676). As illustrated in FIG. 38, the weight computed from the weight function (683 and 685) deceases exponentially as an elapsed time used in the weight function (683 and 685) increases and the weights (e.g., 682 and 684) computed from the weight function (683 and 685) based a first elapsed time (678) before the respective time instance (681) and a second elapsed time (679) after the respective time instance (681) are an equal in magnitude but have opposite signs, when the first elapsed time (678) and the second elapsed time (679) are equal in length.

In one embodiment, to determine whether or not a location having a set of coordinates, in the plurality of locations of the respective mobile device, is located within the predetermined region, the computing device: combines the coordinates into an identifier of a cell among the plurality of cells in the grid reference system, where the cell contains the respective location on the surface of the Earth; and searches a set of cell identifiers pre-associated with the predetermined region to identify a match. The location having the set of coordinates is determined to be within the predetermined region in response to the match being found in the set of cell identifiers pre-associated with the predetermined region; and the location having the set of coordinates is determined to be not within the predetermined region in response to the match not being found in the set of cell identifiers pre-associated with the predetermined region.

For example, the coordinates are combined to provide the identifier of the cell according to a predetermined function of the coordinates of the respective location.

For example, the coordinates are combined by: generating two integers from the coordinates according to a precision level; and combining the two integers to provide the identifier of the cell.

For example, the coordinates are combined by: selecting digits from the coordinates in accordance with a cell resolution level; and combining the digits selected from the coordinates into an integer representing the identifier of the cell.

For example, the method to combine coordinates of a location to a cell identifier as illustrated in FIG. 13 can be used.

The techniques of FIGS. 37-41 are discussed in the application of constructing location responses from visitation measurements. The techniques can be also used in the measurement of other types of responses, such as transaction responses, web visitation responses, service subscription responses, retail store visitation response, etc. For example, in FIG. 37, the location visits (671, 573, . . . , 675) in time segments (e.g., t1, t2, . . . , tn) can be replaced with other types of measurements, such as web visits, payment transactions, purchases, retail store visits, service subscription for the respective time segments in the time interval between time instances Ta and Tb to provide a measurement Vb (677). The anti-symmetric weight function (683 and 685) of FIG. 38 can be applied to such measurements to construct a change response (687) at a time instance (681) as illustrated in FIG. 38; and a series of change responses for a time period (691) can be added up to generate a total response (693) for the time period (691).

Activity Area

After a period of time of influences, a set of mobile devices can be divided into a set of exposed mobile devices that have been exposed to the influence within the period of time and a set of non-exposed mobile devices that have not been exposed to the influence within the period of time. A statistical change in location behavior of mobile devices as a result of the influences can be measured based on differences between the behavior changes of exposed mobile devices and the behavior changes of non-exposed mobile devices, where the behavior changes can be quantified using the methods of FIGS. 37-41.

It is desirable to compare the behavior changes of exposed and non-exposed mobile devices that have similar characteristics or similar distributions of characteristics in the calculation of the statistical change. Similar mobile devices have similar profiles of keywords, which may be determined using the methods of FIGS. 28-36.

For example, a balanced approach can be used in selecting the mobile devices in the statistical calculation of the location behavior such that the population of the exposed mobile devices used in the calculation matches with the population of the non-exposed mobile devices used in the calculation in population size and/or keyword profile distribution.

For example, when the population of the exposed mobile devices is significantly smaller than the population of the non-exposed mobile devices, it is possible to select a subset of the non-exposed mobile devices such that the selected subset matches with the exposed mobile devices in population size and/or keyword profile distribution.

For example, for each of the exposed mobile devices, a corresponding non-exposed mobile device that has a keyword profile best matching the keyword profile of the corresponding exposed mobile device can be selected into the subset used in the statistical calculation. As a result, the selected subset of the non-exposed mobile devices and the set of exposed mobile devices have a matched population size and a matched keyword profile distribution.

Similarly, when the population of the non-exposed mobile devices is significantly smaller than the population of the exposed mobile devices, it is possible to select a subset of the exposed mobile devices such that the subset matches with the set of non-exposed mobile devices in population size and/or keyword profile distribution.

In some instances, when an exposed (or non-exposed) mobile device does not have a corresponding matching non-exposed (or exposed) mobile device according to a predetermined matching threshold, the exposed (or non-exposed) mobile device can be excluded from the statistical calculation. Thus, the matching process may filter both the population of exposed mobile devices and the population of the non-exposed mobile devices to reach a subset of exposed mobile devices and a subset of non-exposed mobile devices, where both subsets have a matching population size and a matching keyword profile distribution. A statistical change in location behavior of mobile devices visiting a region as a result of the time period of influences for the region can be calculated from the differences between the behavior changes of the subset of exposed mobile devices and the behavior changes of the subset of non-exposed mobile devices.

In the balanced approach, the difference between the average responses of exposed mobile devices and the average responses of the non-exposed mobile devices is the same as the average of the response differences between each pair of exposed mobile device and its matching non-exposed mobile device.

However, it is also possible to use an unbalanced approach in selecting the mobile devices in the statistical calculation of the location behavior where the population size of the exposed mobile devices used in the calculation is significantly different from the population size of the non-exposed mobile devices used in the calculation. When an unbalanced approach is used, the statistical change in location behavior of mobile devices as a result of the time period of influences can be measured from the difference in the average response of the exposed mobile devices used in the calculation and the average response of the non-exposed mobile devices used in the calculation.

In one embodiment, the exposed and non-exposed mobile devices used in the statistical calculation are selected according to an activity area relevant to the location behavior. Exposed and non-exposed mobile devices that are not within the activity area are excluded in the calculation or measurement of the statistical change.

FIG. 42 illustrates an activity area for the selection of mobile devices for measuring a statistical change in location behavior of mobile devices according to one embodiment.

In FIG. 42, a region (101) has a predetermined boundary; and a set of influences are provided to exposed mobile devices (e.g., 109, 735) during a period of time. Some of the mobile devices are not exposed to the influences (e.g., 733, 736).

In FIG. 42, the region (101) has an activity area (731). Examples of the activity area (731) include a business area that covers the locations of customers of the region (101), a trade area of the merchant at the region (101). For example, the activity area (731) is such that a customer located outside of the activity area (731) is unlikely to use the service of the region (101). For example, the activity area (731) is such that a customer located outside of the activity area (731) is unlikely to travel to the region (101) to obtain services similar to the service offered by the region (101) (e.g., due to excessive travel time and/or travel distance). For example, the activity area (731) is such that a customer located outside of the activity area (731) is unlikely to travel to a local commercial district (739) (e.g., a strip mall, a shopping center), in which the region (101) is located, to obtain services similar to the services offered by the merchants in the commercial district (739) (e.g., due to excessive travel time and/or travel distance).

In FIG. 42, the activity area (731) is a circular area defined by an activity radius (737) from the location of the region (101) (e.g., the center location, or the location of the main entrance, of the region (101)). The activity radius (737) can be a typical or average radius of business areas of merchants that provide the same or similar service offered by the region (101).

FIG. 42 shows an alternative activity area (732) that has an irregular shape determined by the locations of the customers of the region (101) and/or the local commercial district (739). For example, the activity area (732) can be determined based on identifying a set of mobile devices (e.g., 109, 733) that have visited the region (101) and/or the commercial district (739) during a period of time and the home locations and/or visited locations of the set of mobile devices (e.g., 109, 733); and the alternative activity area (732) is constructed to contain the home locations and/or visited locations (e.g., covering area where the density of such locations is above a threshold). For example, the activity area (732) can be calculated based on the locations from which the mobile devices (e.g., 109, 733) visit the region (101) (e.g., from the home/work locations of the mobile devices (e.g., 109, 733)). For example, the activity area (732) can be determined from a union of cells of a grid reference system (e.g., as illustrated in FIGS. 2-12), where each of the cells have a predetermined resolution/size and have more than a predetermined number of the home locations and/or visited locations.

FIG. 42 illustrates a boundary of a local commercial district (739) using a circle defined by a radius (738). For example, the radius may be the average radius of commercial district (739) in the area, or a predetermined radius (e.g., 1 mile, or half a mile). In other instances, the local commercial district (739) is identified based on the actual boundary of the parcels that are in the vicinity of the region (101) and are allocated for commercial use.

In some instances, the activity area (732) is a combination of a plurality of activity areas (e.g., 731) determined from different approaches discussed above. Different weights can be applied to activity areas (e.g., 731) determined from different approaches; cells of the activity areas (e.g., 731) identified by different approaches receive weight contributions from the respective approaches; and cells having a combined weight above a threshold can be selected to collectively form a resulting combined activity area (732).

In determining the statistical change in location behavior of the mobile devices visiting the region (101), the system uses the location responses of the mobile devices (e.g., 109, 733) that are within the activity area (731) and excludes the location responses of the mobile devices (e.g., 735, 736) that are located outside the activity area (731).

In some instances, when a household or home location of a mobile device (109 or 733) is within the activity area (731), the mobile device (109 or 733) is determined to be within the activity area (731). Some methods to identify the household or home location of a mobile device (109 or 733) are provided in U.S. Pat. No. 9,503,863, entitled “Systems and Methods to Identify Home Addresses of Mobile Devices”, the entire disclosure of which is hereby incorporated herein by reference. Other methods to identify the household or home location of a mobile device (109 or 733) (e.g., based on a registration process and/or a statistical analysis) can also be used.

In other instances or in combination, when a mobile device (109 or 733) has a pattern of visitation locations within the activity area (731) during the period of time of the influence (and/or other time periods), the mobile device (109 or 733) is determined to be within the activity area (731).

FIG. 43 shows a method to measure a statistical change in location behavior of mobile devices based on an activity area according to one embodiment. For example, the method of FIG. 43 can be used to select the exposed and non-exposed mobile devices according to the activity area illustrated in FIG. 42 in computing the statistical measurement of changes in FIG. 41.

In FIG. 43, a computing device (e.g., server (187) illustrated in FIG. 1) is configured to: receive (741) an identification of a region (101) having a predetermined boundary; identify (743) an activity area (731 or 732) containing the region (101); receive (745) exposure statuses of tracked mobile devices (e.g., 109, 733, 735, . . . , 736) with respect to a time period of influences to visit the region (101); identify (747) a first subset of the tracked mobile devices that reside within the activity area (731 or 732) and are exposed to the influences (according to the exposure statuses) (e.g., 109, 735); identify (749) a second subset of the tracked mobile devices that reside within the activity area (731 or 732) but are not exposed to the influences (e.g., 733, 736); compute (751), for each mobile device in the first and second subsets (e.g., 109, 733), a change of visitation to the region (101) with respect to the time period; and compute (753) a statistical measurement of changes resulting from the time period of influences from a difference between the average of the visitation changes of mobile devices (e.g., 109) in the first subset and the average of the visitation changes of mobile devices in the second subset (e.g., 733).

For example, for each of the mobile devices (e.g., 109, 733) in the first and second subsets, a change response (693) illustrated in FIG. 39 is computed using the method of FIG. 40; and a statistical measurement of changes resulting from the time period of influences is computed from a difference between the average response of the first subset of exposed mobile devices (e.g., 109) and the average response of the second subset of non-exposed mobile devices (e.g. 733) in a way similar to the method of FIG. 41.

In general, the first subset of exposed mobile devices and the second subset of non-exposed mobile devices have different population sizes. However, since the first subset of exposed mobile devices and the second subset of non-exposed mobile devices are selected based on the activity area (731 or 732), the first subset of exposed mobile devices and the second subset of non-exposed mobile devices are likely to have similar keyword profiles in an average sense, especially when the population sizes of the first and second subsets are sufficiently large.

Alternatively, a balance approach can be applied to the first subset of exposed mobile devices and the second subset of non-exposed mobile devices so that the exposed mobile devices selected from the first subset and the non-exposed mobile devices selected from the second subset for the computation of the statistical measurement have the same population size and matched keyword profiles.

In FIG. 43, the activity area (731 or 732) can be determined using one of the methods discussed in connection with FIG. 42 or one of the methods of FIG. 44-46.

FIG. 44 shows a method to identify an activity area. In FIG. 44, the method includes: receiving (761) an identification of a region (101) having a predetermined boundary; identifying (763) a representative location of the region (101); and determining (765), based on a predetermined radius (737) from the representative location, an activity area (731) for the selection of mobile devices (e.g., 109, 733) in determining changes of visitation to the region (101).

In some instances, the radius (737) can be replaced with a threshold travel distance on a road system from the representative location of the region (101) (e.g., the center of the region (101), an entrance of the region (101), a street address of the region (101)).

In some instances, the radius (737) can be replaced with a threshold travel time on a road system from the representative location of the region (101) (e.g., the center of the region (101), an entrance of the region (101), a street address of the region (101)) (e.g., based on the travel speed(s) according to the speed limit(s) of the roads).

FIG. 45 shows another method to identify an activity area. In FIG. 45, the method includes: receiving (771) an identification of a region (101) of a merchant location having a predetermined boundary; identifying (773) transactions of users with the merchant (or the merchants in the local commercial district (739)); and determining, (775) based on household locations of the users, an activity area (732) for the selection of mobile devices in determining changes of visitation to the region (101).

FIG. 46 shows a further method to identify an activity area. In FIG. 46, the method includes: receiving (781) an identification of a first region (101) of a merchant location having a predetermined boundary; identifying (783) a plurality of merchants that are located within a second region (739) and are adjacent to the merchant at the first region (101); identifying (785) mobile devices that have been to the second region (739) and/or having transactions (or interactions) with the plurality of merchants within a period of time (which may or may not be the same as, and may or may not have overlap with, the time period of influence); determining (787) household locations (and/or visited locations) of the mobile devices within the period of time; and determining (789), based on household locations (and/or visited locations) of the mobile devices, an activity area (732) for the selection of mobile devices in determining changes of visitation to the first region (101). The activity area (732) is determined to cover regions where the household locations (and/or visited locations) of the mobile devices have a density above a threshold.

The methods of FIGS. 44-46 can be combined via a weighted average method. For example, the activity areas (e.g., 731, 732) identified using the methods of FIGS. 44-46 (and/or other methods) can be mapped to corresponding set of cells of a grid reference system, where the cells are contained within the respective activity areas. Each of the methods used to identify an activity area is assigned a weight; and when a cell is within the activity area identified by a method, the weight of the method is added to the cell. Cells that have an accumulated weight above a threshold are selected as part of the combined activity area generated from the weighted average method.

For example, a system to measure a statistical change in location behavior of mobile devices may include a plurality of mobile devices (e.g., 109, 733, 735, 736) and a computing device (e.g., server (187)) in communication with the mobile devices (e.g., 109, 733, 735, 736) to receive locations of the mobile devices (e.g., 109, 733, 735, 736) at different time instances. Each mobile device (e.g., 109) of the plurality of mobile devices (e.g., 109, 733, 735, 736) determines locations of the mobile device (e.g., 109) at respective time instances when the mobile device (e.g., 109) is at the respective locations (e.g., 111).

After receiving information identifying a time period (e.g., a day, a week, or a month) during which at least one message related to a predetermined region (101) is transmitted to a portion of the mobile devices (e.g., 109, 733, 735, 736) at least once, the computing device (e.g., server (187)) determines an activity area (731 or 732) associated with the predetermined region (101) (e.g., using one of the methods of FIGS. 44-46).

Based on the activity area (731 or 732), the computing device (e.g., server (187)) identifies: a first subset of exposed mobile devices (e.g., 109) that are in the activity area (731 or 732) during the time period and that have received the message at least once during the time period; and a second subset of non-exposed mobile devices (e.g., 733) that are also in the activity area (731 or 732) during the time period but have not received the message during the time period.

Typically, the first subset and the second subset have different population sizes; and the first subset and the second subset are selected without individually matching mobile devices in the first subset to mobile devices in the second subset based on keyword profiles of mobile devices. The plurality of mobile devices (e.g., 109, 733, 735, 736) generally include at least one mobile device (e.g., 735) that is not in the activity area during the time period and receives the message at least once during the time period and at least one mobile device (e.g., 736) that is not in the activity area during the time period and does not receive the message during the time period.

Optionally, the computing device determines whether a particular mobile device (e.g., 109) is within the activity area (731 or 732) during the time period based on whether the particular mobile device has a home location or a location visitation pattern in the activity area (731 or 732) during the time period.

For each respective mobile device (e.g., 109 or 733) in the first subset and the second subset, the computing device (e.g., server (187)) determines a response (e.g., 693) of the respective mobile device (e.g., 109 or 733) based on counts of visitation to the predetermined region (101) by the respective mobile device (e.g., 109 or 733).

The computing device (e.g., server (187) then calculates a measurement of a location behavior change based on a difference between: an average of responses of mobile devices in the first subset (e.g., 109), and an average of responses of mobile devices in the second subset (e.g., 733).

The activity area (e.g., 731 or 732) can be determined to cover locations of first mobile devices that have visited the predetermined region (101) within a second period of time (e.g., a week, a month, three months, half a year), such as the home locations of the first mobile devices during the second period of time and/or the visited locations of the first mobile devices during the second period of time.

Alternatively, or in combination, the activity area (e.g., 731 or 732) can be determined to cover the locations of the first mobile devices that have visited a surrounding area (739) of the predetermined region (101) within the second period of time. The second period of time may or may not include the time period (e.g., a day, a week, or a month) during which the message related to the predetermined region (101) (101) is transmitted to mobile devices (e.g., 109, 735).

In some instances, the activity area (e.g., 731 or 732) is determined based on a representative location of the predetermined region (101) and a threshold radius, a threshold travel distance from the representative location on a road system, and/or a threshold travel time from the representative location on the road system.

Optionally, the counts of visitation to the predetermined region (101) by a particular mobile device (109) are determined based on whether or not each location of the mobile device (109), as represented by a set of coordinates, is located within the predefined boundary of the predetermined region (101). Such a determination can be performed by the computing device (e.g., server (187)) by: combining the coordinates into an identifier of a cell among the plurality of cells in the grid reference system; and searching a set of cell identifiers pre-associated with the predetermined region to identify a match, as illustrated in FIG. 12. The location having the set of coordinates is determined to be within the predetermined region in response to the match being found in the set of cell identifiers pre-associated with the predetermined region; and the location having the set of coordinates is determined to be not within the predetermined region in response to the match not being found in the set of cell identifiers pre-associated with the predetermined region. The coordinates can be combined to provide the identifier of the cell according to a predetermined function of the coordinates of the respective location, as illustrated in FIG. 12 or FIG. 13.

Matching

One embodiment disclosed in the present application matches exposed mobile devices and non-exposed mobile devices having similar behaviors in an attribute space to measure the effect of the exposure based on clusters of mobile devices in the attributes.

In general, the distribution of exposed mobile devices and the distribution of non-exposed mobile devices may not agree with each other in a given population of exposed and non-exposed mobile devices. Thus, when the average visitation change response of the exposed mobile devices in the population is compared directly to the average visitation change response of the non-exposed mobile devices in the population, the difference may not accurately reflect the comparison between similar mobile devices.

Alternatively, a selection process can be used to individually paring an exposed mobile device with its most similar non-exposed mobile device in an attribute space. However, paring individual mobile devices is computationally inefficient.

When paring exposed mobile devices and non-exposed mobile device based on clusters in attribute space, the distribution of the exposed mobile devices matches with the distribution of non-exposed mobile devices in a statistical way. Thus, the combined performance of computation efficiency and statistical accuracy of the exposure influence measurement is improved.

FIG. 47 illustrates clustering mobile devices in an attribute space to match exposed mobile devices and non-exposed mobile devices.

In FIG. 47, each dot (e.g., 802 or 804) represents a mobile device (e.g., 109) in an attribute space. A filled dot (e.g., 802) represents a mobile device having one exposure status (e.g., exposed); and a non-filled dot (e.g., 804) represents a mobile device having another exposure status (e.g., non-exposed).

Attributes used to characterize the mobile devices in the attribute space may include keyword values computed using the techniques of FIGS. 28-36, the relative position of the mobile devices in relation with the activity area of a region as discussed in FIGS. 42-46, and/or other attributes of the mobile devices, such as demographic attributes, behavioral attributes, propensity attributes, and/or psycho-graphic attributes that are computed or received using various methods.

The values of the attributes for a given mobile device characterize the behavior of the given mobile device. Thus, two mobile devices that are close to each other in the attribute space are similar in their behaviors.

Through a cluster analysis, mobile devices that are close to each other in the attribute space can be identified as a cluster (e.g., 801, 803, 805, or 807).

For example, the cluster analysis can be performed using any of the k-mean clustering algorithms known in the field; and each cluster can be labeled using a unique reference number that can also be considered as a score of the mobile devices in the cluster.

The clusters can also be identified by scoring the mobile devices according to their attribute values (e.g., using a continuous function that maps the attribute values to a score). Mobile devices having the same score, or scores in the same range, can be considered as in the same cluster.

As illustrated in FIG. 47, some of the clusters (e.g., 803 or 805) many be homogenous in status for containing only mobile devices that have the same exposure status. For example, one cluster (803) has only exposed mobile devices; and another cluster (805) has only non-exposed mobile devices. Thus, meaningful comparisons cannot be made for the mobile devices in such clusters (e.g., 803 and 805); and such clusters are eliminated from the computation of the measurement of the influence of the exposure.

After eliminating the clusters that contain only one kind of mobile devices (either exposed or non-exposed), the remaining mobile devices in the population is illustrated in FIG. 48.

FIG. 48 illustrates a selection of unbalanced matching of exposed and non-exposed mobile devices in clusters. For example, the clusters (801 and 807) illustrated in FIG. 48 can be obtained by removing clusters that only contains mobile devices having the same exposure status (e.g., either exposed or non-exposed) and/or removing clusters having population sizes for a given exposure status that are below a threshold value.

In FIG. 48, each cluster may have different numbers of mobile devices that have different exposure statuses. For example, in one cluster (e.g., 801), the population size of non-exposed mobile devices (e.g., 804) in the cluster (801) is larger than the population size of exposed mobile devices (e.g., 802) in the cluster (801); and in another cluster (e.g., 807), the population size of non-exposed mobile devices (e.g., 808) in the cluster (807) is smaller than the population size of exposed mobile devices (e.g., 806) in the cluster (807). Thus, the populations of exposed and non-exposed mobile devices in at least some clusters do not match with each other in size.

When the populations of exposed and non-exposed mobile devices in individual clusters do not match with each other in size, the effect of exposure can be computed based on the cluster-level difference between the cluster-level average response of the exposed mobile devices and the cluster-level average response of the non-exposed mobile devices in each cluster (e.g., 801 and 807). Cluster-level differences can be combined to generate a measurement for the set of clusters (e.g., 801 and 807).

For example, for each cluster (e.g., 801 or 807), the responses (e.g., 693) of exposed mobile devices (e.g., 802 or 806) computed based on counts of visitation to the predetermined region (101) by the respective mobile devices (e.g., 802 or 806) (e.g., using the technique of FIGS. 37-39) can be added up and then divided by the number of exposed mobile devices (e.g., 802 or 806) in the cluster (e.g., 801 or 807) to compute the average response of exposed mobile devices (e.g., 802 or 806) of the cluster (e.g., 801 or 807). Similarly, for each cluster (e.g., 801 or 807), the responses (e.g., 693) of non-exposed mobile devices (e.g., 804 or 808) computed based on counts of visitation to the predetermined region (101) by the respective mobile devices (e.g., 804 or 808) can be added up and then divided by the number of non-exposed mobile devices (e.g., 804 or 808) in the cluster (e.g., 801 or 807) to compute the average response of non-exposed mobile devices (e.g., 804 or 808) of the cluster (e.g., 801 or 807). The difference between the average response of exposed mobile devices (e.g., 802 or 806) of the cluster (e.g., 801 or 807) and average response of non-exposed mobile devices (e.g., 804 or 808) of the cluster (e.g., 801 or 807) provides a measurement of the influence of the exposure for the cluster (e.g., 801 or 807).

The influence measurements of the exposure computed for different clusters (e.g., 801, 807) can be summed up for an average measurement of the influence of the exposure on mobile devices in general. Preferably, the average is computed based on weights corresponding to the sizes of the clusters. For example, the influence measurement of each cluster (e.g., 801) is multiplied by a size measurement of the cluster (e.g., 801) to obtain the size-weighted influence measurement; and the size-weighted influence measurements for the clusters (e.g., 801, 807) are summed up and divided by the sum of the size measurements of the clusters (e.g., 801, 807) to obtain a weighted average of influence measurements. A size measurement of a cluster (801) can be, for example, proportional to a count of exposed mobile devices in the cluster (e.g., 801), a count of non-exposed mobile devices in the cluster (e.g., 801), or a count of all mobile devices, regardless of their exposure statuses, in the cluster (e.g., 801).

In some instances, the system eliminates not only clusters that have no exposed mobile device and clusters that have no non-exposed mobile device, but also clusters each having less than a threshold count of exposed mobile devices and clusters each having less than a threshold count of non-exposed mobile devices. Thus, the average responses computed for the remaining clusters are statistically significant.

In some instances, the system down-samples the clusters to match the populations of exposed mobile devices and non-exposed mobile devices in each cluster, as illustrated in FIG. 49.

FIG. 49 illustrates a selection of balanced matching of exposed and non-exposed mobile devices in clusters. For example, the clusters (801′ and 807′) illustrated in FIG. 49 can be obtained by down-sampling the clusters (801 and 807) illustrated in FIG. 48 to match the population sizes of exposed mobile devices and non-exposed mobile devices in each cluster.

For example, the non-exposed population (e.g., mobile device (804)) in a cluster (801) is larger than the exposed population (e.g., mobile device (802)) in the cluster (801), a subset of non-exposed population can be selected (e.g., randomly) and removed from the cluster (801) to generate the cluster (801′) that has matched population sizes between exposed mobile devices (e.g., 802) and non-exposed mobile devices (e.g., 804). Thus, the same number of non-exposed mobile devices (e.g., 804) are selected (e.g., randomly) from the cluster (801) to match the number of exposed mobile devices (e.g., 802) in the cluster (801) to generate the down-sized cluster (801′).

For example, the exposed population (e.g., mobile device (806)) in a cluster (807) is larger than the non-exposed population (e.g., mobile device (808)) in the cluster (807), a subset of exposed population can be selected (e.g., randomly) and removed from the cluster (807) to generate the cluster (807′) that has matched population sizes between exposed mobile devices (e.g., 806) and non-exposed mobile devices (e.g., 808). Thus, the same number of non-exposed mobile devices (e.g., 808) are selected (e.g., randomly) from the cluster (807) to match the number of exposed mobile devices (e.g., 806) in the cluster (807) to generate the down-sized cluster (807′).

When each of the down-sized clusters (801′ and 807′) has a matched population size between exposed and non-exposed mobile devices, as illustrated in FIG. 49, it is not necessary to compute the differences of average responses of exposed and non-exposed mobile devices at the cluster level and then perform size-weighted average across the clusters (801′ and 807′). Weighting according to the cluster sizes measured by counting the exposed mobile devices, non-exposed mobile devices, or both, provides the same result as taking the average of mobile devices across the clusters.

FIG. 50 shows a method to compute a visitation change measurement using balanced matching of mobile devices. For example, the method of FIG. 50 can be implemented in the system illustrated in FIG. 1.

In FIG. 50, a computing device (e.g., the server (187)) is configured to perform (811) a cluster analysis of a population of mobile devices (e.g., as illustrated in FIG. 47) in an attribute space, where each of the mobile devices has either a first status or a second status (e.g., exposed or non-exposed). The cluster analysis identifies a plurality of clusters of mobile devices in the attribute space.

For example, the attribute space may include attributes corresponding to the keyword values computed using the techniques of FIGS. 28-36. Further, the attribute space may further include an attribute identifying whether or a mobile device is within an activity area discussed in FIGS. 42-46.

After the cluster analysis (811), the computing device is configured to: identify (813), among the plurality of clusters (e.g., 801, 803, 805, 807), a cluster (e.g., 801 or 807) that has at least one mobile device having the first status and at least one mobile device having the second status; and select, (815) from mobile devices in the cluster, a number of first mobile devices having the first status (e.g., exposed) and the same number of second mobile devices having the second status (e.g., non-exposed).

For example, if the cluster (807) has more exposed mobile devices than the non-exposed mobile devices, the computing device may select all of the non-exposed mobile devices and a number of exposed mobile devices that has a population size equal to the population size of the non-exposed mobile devices in the cluster (807). Alternatively, the computing device may select some of the non-exposed mobile devices and a number of exposed mobile devices that has a population size equal to the population size of the selected non-exposed mobile devices in the cluster (807).

If it is determined (817) that there is another cluster from which mobile devices can be selected, the operations (813 and 815) are repeated for the next cluster.

After all of the clusters are processed and an equal number of first mobile devices having the first status (e.g., exposed) and second mobile devices having the second status (e.g., non-exposed) are selected, the computing device is configured to compute (819) a visitation change measurement from an average of visitation changes of a number of the selected first mobile devices having the first status (e.g., exposed) in the identified clusters and an average of visitation changes of the same number of the selected second mobile devices having the second status (e.g., non-exposed) in the identified clusters.

For example, the visitation change measurement can be computed as the difference between the average of visitation changes of the number of the selected first mobile devices having the first status (e.g., exposed) in the identified clusters (e.g., 801′ and 807′) and the average of visitation changes of the same number of the selected second mobile devices having the second status (e.g., non-exposed) in the identified clusters (e.g., 801′ and 807′), where the visitation change of each mobile device is in the form of a response (693) computed using the technique of FIGS. 37-39.

Down-sampling the populations for a match between mobile devices of different statuses (e.g., exposed and non-exposed) may significantly alter the distribution of the mobile devices. As a result, the visitation change measurement, computed from the down-sampled populations to create the perfect match between the exposed and non-exposed, may not accurately reflect the population distribution in general.

In some instances, the distribution of the down-sampled population of the mobile devices is compared to the distribution of the entire population of the mobile devices for the assessment of the quality of the down-sampling operation. When the distribution of the down-sampled population of the mobile devices deviates significantly from the distribution of the entire population of the mobile devices, it may be preferred to compute the visitation change measurement without requiring the perfect population match at cluster level. For example, the deviation in distribution can be qualified and examined by computing and comparing the entropy of the distributions of the mobile devices in the attribute space.

In some instances, when the exposed population is larger than the non-exposed population in a cluster (e.g., 807), the exposed population is not down-sampled to use an unbalanced cluster (e.g., 807). In other instances, both the exposed population and the non-exposed population in a cluster (e.g., 801 or 807) are not down-sampled, as illustrated in FIG. 48. The method of FIG. 51 can be used to compute the visitation change measurement for clusters (e.g., 801 and 807) that have unmatched population sizes for exposed and non-exposed.

FIG. 51 shows a method to compute a visitation change measurement using balanced or unbalanced matching of mobile devices. For example, the method of FIG. 51 can be implemented in the system illustrated in FIG. 1.

In FIG. 51, a computing device (e.g., the server (187)) is configured to perform (821) a cluster analysis of a population of mobile devices (e.g., as illustrated in FIG. 47), where each of the mobile devices has either a first status or a second status regarding to exposure to period of influence related to a region (101) (e.g., exposed or non-exposed). The cluster analysis identifies a plurality of clusters of mobile devices in an attribute space, in a way similar to the cluster analysis (811) of FIG. 50.

For example, the attribute space may include attributes corresponding to the keyword values computed using the techniques of FIGS. 28-36. Further, the attribute space may optionally include an attribute identifying whether or a mobile device is within an activity area discussed in FIGS. 42-46.

After the cluster analysis (821), the computing device is configured to: identify (823), among the plurality of clusters (e.g., 801, 803, 805, 807), a cluster (801 or 807) that has at least one mobile device having the first status and at least one mobile device having the second status; and compute (825), for the cluster, a difference between an average of changes in visitation to the region (101) by first mobile device having the first status in the cluster (801 or 807) and an average of changes in visitation to the region (101) by second mobile devices having the second status in the cluster (801 or 807).

In general, the number of the first mobile device having the first status (e.g., exposed) in the cluster (801 or 807) is not equal to the number of the second mobile devices having the second status (e.g., non-exposed) in the cluster (801 or 807). However, in some instances, the number of the first mobile device having the first status in the cluster (801′ or 807′) is the same as the number of the second mobile devices having the second status in the cluster (801′ or 807′) (e.g., as a result of an optional down-sampling operation, or by chance).

If it is determined (827) that there is another cluster that has mobile devices having different exposure status and thus can be selected to compute (825) a difference at the cluster level, the operations (823 and 825) are repeated for the next cluster.

After all of the clusters are processed, the computing device computes (829) a visitation change measurement from an average of the differences computed for the clusters. For example, the average can be computed based on weights proportional to the sizes of the clusters (e.g., the population size of the exposed mobile devices, the population size of the non-exposed mobile devices, or the population size of both exposed and non-exposed mobile devices in a corresponding cluster).

FIG. 52 shows a method to compute a visitation change measurement based on matching population distributions. For example, the method of FIG. 52 can be implemented in the system illustrated in FIG. 1, in accordance with the method of FIG. 50 or 51.

In FIG. 52, a computing device (e.g., the server (187) or another computing apparatus coupled to the server (187)) is configured via instructions to: receive (831) location information of a plurality of mobile devices (e.g., 109), each determining its locations at respective time instances; receive (833) information identifying a time period (691) during which at least one message related to a predetermined region (101) is transmitted to a portion of the mobile devices (e.g., 109) at least once; identify (835) a plurality of clusters (801, 803, 805, 807) of the plurality of mobile devices in an attribute space (e.g., as illustrated in FIG. 47); select (837) a subset of the plurality of mobile devices based at least in part on identifications of the clusters in the attribute space (e.g., a subset of devices illustrated in FIG. 48 or 49); identify (839), in the subset, a plurality of first mobile devices each having a first status with respect to the message (e.g., exposed mobile devices (e.g., 802 and 806)); identify, (841) in the subset, a plurality of second mobile devices each having a second status with respect to the message (e.g., non-exposed mobile devices (e.g., 804 and 808); for each respective mobile device (e.g., 802, 804, 806, or 808) in the first mobile devices and the second mobile devices, compute (843) a response (693) of the respective mobile device based on counts of visitation to the predetermined region by the respective mobile device (e.g., as illustrated in FIGS. 37-39); and calculate (845) a measurement of a location behavior change based on a difference between responses of the first mobile devices and responses of the second mobile devices.

For example, each response of a mobile device can be computed based on applying, at a time instance (e.g., 681), an anti-symmetric weight function (683 and 685) having an exponential distribution to a time series of visitation measurements (e.g., 677, 676 illustrated in FIG. 38) to generate a response (687) corresponding to the time instance (681), and summing the responses (687, 688, . . . , 689) corresponding to time instances falling within the predetermined period of time (691) as the response (693 illustrated in FIG. 39) of the mobile device.

For example, for each respective cluster (e.g., 801) that contains at least a portion of the first mobile devices and the second mobile devices, a cluster-level measurement of the location behavior change is computed from a difference between: a first average of responses of devices that are both in the plurality of the first mobile devices having the first status (e.g., exposed) in the selected subset and in the respective cluster (e.g., 801); and a second average of responses of devices that are both in the plurality of the second mobile devices having the second status (e.g., non-exposed) in the selected subset and in the respective cluster (e.g., 801). The cluster-level measurements of the location behavior change of clusters (e.g., 801 and 807) that contain mobile devices in the subset are combined to generate the measurement of the location behavior change.

For example, the cluster-level measurements can be combined based on a weighted average where each cluster-level measurement for a cluster is given a weight proportional to the size of the cluster. The size of the cluster can be computed based on the number of mobile devices having the first status (e.g., exposed) in the cluster, the number of mobile devices having the second status (e.g., non-exposed) in the cluster, or the number of mobile devices selected in the subset regardless of their statuses.

For example, the first mobile devices can be identified to have the first status (e.g., exposed) for having received the message during the time period; and the second mobile devices can be identified to have the second status (e.g., non-exposed) for having not received the message during the time period. The subset of the plurality of mobile devices is selected by at least excluding mobile devices in one or more clusters that are homogeneous in status (e.g., a cluster having either mobile devices all having the first status or mobile devices all having the second status). In some instances, clusters each having mobile devices of a same status (e.g., exposed or non-exposed) less than a threshold population size are also excluded from the subset.

In general, the population size of the first mobile devices in the subset having the first status (e.g., exposed) is different from the population size of the second mobile devices in the subset having the first status (e.g., non-exposed). The unmatched population sizes can also occur in one or more clusters of the selected subset of mobile devices, as illustrated in FIG. 48.

Optionally, random down-sampling can be performed at cluster level to balance the population sizes of the selected mobile devices in the subset in each cluster, such that for each cluster of the subset of mobile devices, the population size of mobile devices having the first status (e.g., exposed) matches with the population size of mobile devices having the second status (e.g., non-exposed), as illustrated in FIG. 49.

When the population sizes are balanced at cluster level, it is not necessary to first compute the measurements at the cluster-level and then combine the cluster-level measurements; and the average of the mobile devices in the selected subset and having the same status can be computed and the averages for different statuses be compared to compute the measurement for the selected subset. However, the balanced approach requires a down-sampling operation, which may further change the distribution of the selected mobile devices for the computation and thus less accurate in reflecting the measurement for the population that is not down-sampled.

In computing (843) of the response of the respective mobile device, the counts of visitation to the predetermined region by the respective mobile device can be determined based on whether or not each location, in a plurality of locations of the respective mobile device at a plurality of time instances respectively, is located within the predetermined region.

To determine whether or not a location having a set of coordinates, in the plurality of locations of the respective mobile device, is located within the predetermined region, the techniques of FIGS. 1-14 can be used. For example, the computing device may combine the coordinates into an identifier of a cell among the plurality of cells in the grid reference system, where the cell contains the respective location on the surface of the Earth; and search a set of cell identifiers pre-associated with the predetermined region to identify a match. The location having the set of coordinates is determined to be within the predetermined region in response to the match being found in the set of cell identifiers pre-associated with the predetermined region; and the location having the set of coordinates is determined to be not within the predetermined region in response to the match not being found in the set of cell identifiers pre-associated with the predetermined region.

For example, the coordinates can be combined to provide the identifier of the cell according to a predetermined function of the coordinates of the respective location.

For example, the coordinates can be combined by: generating two integers from the coordinates according to a precision level; and combining the two integers to provide the identifier of the cell, as illustrated in FIG. 12.

For example, the coordinates can be combined by: selecting digits from the coordinates in accordance with a cell resolution level; and combining the digits selected from the coordinates into an integer representing the identifier of the cell, as illustrated in FIG. 13.

Matching the exposed mobile devices and non-exposed mobile devices at cluster-level is computationally more efficient than matching exposed mobile devices and non-exposed mobile devices individually. Further, it allows the use of an unbalance approach, where exposed mobile devices and non-exposed mobile devices participate in the computation do not have to match with each other in population size, as a whole and/or at the cluster level.

The present disclosure includes above disclosed methods implemented in computing devices, computing devices implemented herein the disclosed methods, and the non-transitory computer storage media storing instructions which when executed by computing devices, cause the computing devices to perform the methods. The computing devices can be implemented using one or more data processing systems illustrated in FIG. 15.

FIG. 15 illustrates a data processing system according to one embodiment. While FIG. 15 illustrates various components of a computer system, it is not intended to represent any particular architecture or manner of interconnecting the components. One embodiment may use other systems that have fewer or more components than those shown in FIG. 15.

In FIG. 15, the data processing system (200) includes an inter-connect (201) (e.g., bus and system core logic), which interconnects one or more microprocessors (203) and memory (204). The microprocessor (203) is coupled to cache memory (209) in the example of FIG. 15.

In one embodiment, the inter-connect (201) interconnects the microprocessor(s) (203) and the memory (204) together and also interconnects them to input/output (I/O) device(s) (205) via I/O controller(s) (207). I/O devices (205) may include a display device and/or peripheral devices, such as mice, keyboards, modems, network interfaces, printers, scanners, video cameras and other devices known in the art. In one embodiment, when the data processing system is a server system, some of the I/O devices (205), such as touch screens, printers, scanners, mice, and/or keyboards, are optional.

In one embodiment, the inter-connect (201) includes one or more buses connected to one another through various bridges, controllers and/or adapters. In one embodiment the I/O controllers (207) include a USB (Universal Serial Bus) adapter for controlling USB peripherals, and/or an IEEE-1394 bus adapter for controlling IEEE-1394 peripherals.

In one embodiment, the memory (204) includes one or more of: ROM (Read Only Memory), volatile RAM (Random Access Memory), and non-volatile memory, such as hard drive, flash memory, etc.

Volatile RAM is typically implemented as dynamic RAM (DRAM) which requires power continually in order to refresh or maintain the data in the memory. Non-volatile memory is typically a magnetic hard drive, a magnetic optical drive, an optical drive (e.g., a DVD RAM), or other type of memory system which maintains data even after power is removed from the system. The non-volatile memory may also be a random access memory.

The non-volatile memory can be a local device coupled directly to the rest of the components in the data processing system. A non-volatile memory that is remote from the system, such as a network storage device coupled to the data processing system through a network interface such as a modem or Ethernet interface, can also be used.

In this description, some functions and operations are described as being performed by or caused by software code to simplify description. However, such expressions are also used to specify that the functions result from execution of the code/instructions by a processor, such as a microprocessor.

Alternatively, or in combination, the functions and operations as described here can be implemented using special purpose circuitry, with or without software instructions, such as using Application-Specific Integrated Circuit (ASIC) or Field-Programmable Gate Array (FPGA). Embodiments can be implemented using hardwired circuitry without software instructions, or in combination with software instructions. Thus, the techniques are limited neither to any specific combination of hardware circuitry and software, nor to any particular source for the instructions executed by the data processing system.

While one embodiment can be implemented in fully functioning computers and computer systems, various embodiments are capable of being distributed as a computing product in a variety of forms and are capable of being applied regardless of the particular type of machine or computer-readable media used to actually effect the distribution.

At least some aspects disclosed can be embodied, at least in part, in software. That is, the techniques may be carried out in a computer system or other data processing system in response to its processor, such as a microprocessor, executing sequences of instructions contained in a memory, such as ROM, volatile RAM, non-volatile memory, cache or a remote storage device.

Routines executed to implement the embodiments may be implemented as part of an operating system or a specific application, component, program, object, module or sequence of instructions referred to as “computer programs.” The computer programs typically include one or more instructions set at various times in various memory and storage devices in a computer, and that, when read and executed by one or more processors in a computer, cause the computer to perform operations necessary to execute elements involving the various aspects.

A machine readable medium can be used to store software and data which when executed by a data processing system causes the system to perform various methods. The executable software and data may be stored in various places including for example ROM, volatile RAM, non-volatile memory and/or cache. Portions of this software and/or data may be stored in any one of these storage devices. Further, the data and instructions can be obtained from centralized servers or peer to peer networks. Different portions of the data and instructions can be obtained from different centralized servers and/or peer to peer networks at different times and in different communication sessions or in a same communication session. The data and instructions can be obtained in entirety prior to the execution of the applications. Alternatively, portions of the data and instructions can be obtained dynamically, just in time, when needed for execution. Thus, it is not required that the data and instructions be on a machine readable medium in entirety at a particular instance of time.

Examples of computer-readable media include but are not limited to recordable and non-recordable type media such as volatile and non-volatile memory devices, read only memory (ROM), random access memory (RAM), flash memory devices, floppy and other removable disks, magnetic disk storage media, optical storage media (e.g., Compact Disk Read-Only Memory (CD ROM), Digital Versatile Disks (DVDs), etc.), among others. The computer-readable media may store the instructions.

The instructions may also be embodied in digital and analog communication links for electrical, optical, acoustical or other forms of propagated signals, such as carrier waves, infrared signals, digital signals, etc. However, propagated signals, such as carrier waves, infrared signals, digital signals, etc. are not tangible machine readable medium and are not configured to store instructions.

In general, a machine readable medium includes any mechanism that provides (i.e., stores and/or transmits) information in a form accessible by a machine (e.g., a computer, network device, personal digital assistant, manufacturing tool, any device with a set of one or more processors, etc.).

In various embodiments, hardwired circuitry may be used in combination with software instructions to implement the techniques. Thus, the techniques are neither limited to any specific combination of hardware circuitry and software nor to any particular source for the instructions executed by the data processing system.

The description and drawings are illustrative and are not to be construed as limiting. Numerous specific details are described to provide a thorough understanding. However, in certain instances, well known or conventional details are not described in order to avoid obscuring the description. References to one or an embodiment in the present disclosure are not necessarily references to the same embodiment; and, such references mean at least one.

The use of headings herein is merely provided for ease of reference, and shall not be interpreted in any way to limit this disclosure or the following claims.

Reference to “one embodiment” or “an embodiment” means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment of the disclosure. The appearances of the phrase “in one embodiment” in various places in the specification are not necessarily all referring to the same embodiment, and are not necessarily all referring to separate or alternative embodiments mutually exclusive of other embodiments. Moreover, various features are described which may be exhibited by one embodiment and not by others. Similarly, various requirements are described which may be requirements for one embodiment but not other embodiments. Unless excluded by explicit description and/or apparent incompatibility, any combination of various features described in this description is also included here. For example, the features described above in connection with “in one embodiment” or “in some embodiments” can be all optionally included in one implementation, except where the dependency of certain features on other features, as apparent from the description, may limit the options of excluding selected features from the implementation, and incompatibility of certain features with other features, as apparent from the description, may limit the options of including selected features together in the implementation.

In the foregoing specification, the disclosure has been described with reference to specific exemplary embodiments thereof. It will be evident that various modifications may be made thereto without departing from the broader spirit and scope as set forth in the following claims. The specification and drawings are, accordingly, to be regarded in an illustrative sense rather than a restrictive sense. 

What is claimed is:
 1. A system, comprising: a plurality of mobile devices, wherein each mobile device in the plurality of mobile devices determines locations of the mobile device at respective time instances when the mobile device is at the respective locations; and a computing device in communication with the mobile devices to receive locations of the mobile devices at different time instances, wherein the computing device: receives information identifying a time period during which at least one message related to a predetermined region is transmitted to a portion of the mobile devices at least once; identifies a plurality of clusters of the plurality of mobile devices according to a set of attributes, wherein mobile devices in each of plurality of clusters have values of the attributes that are close to each other in an attribute space defined by the set of attributes; selects a subset of the plurality of mobile devices based at least in part on identifications of the clusters in the attribute space; identifies a plurality of first mobile devices in the subset, wherein each mobile device of the first mobile devices having a first status with respect to the message; identifies a plurality of second mobile devices in the subset, wherein each mobile device of the second mobile devices having a second status with respect to the message; for each respective mobile device in the first mobile devices and the second mobile devices, computes a response of the respective mobile device based on counts of visitation to the predetermined region by the respective mobile device; and calculates a measurement of a location behavior change based on a difference between responses of the first mobile devices and responses of the second mobile devices.
 2. The system of claim 1, wherein the counts of visitation to the predetermined region by the respective mobile device are determined based on whether or not each location, in a plurality of locations of the respective mobile device at a plurality of time instances respectively, is located within the predetermined region.
 3. The system of claim 2, wherein to determine whether or not a location having a set of coordinates, in the plurality of locations of the respective mobile device, is located within the predetermined region, the computing device: combines the coordinates into an identifier of a cell among the plurality of cells in the grid reference system, wherein the cell contains the respective location on the surface of the Earth; and searches a set of cell identifiers pre-associated with the predetermined region to identify a match; wherein the location having the set of coordinates is determined to be within the predetermined region in response to the match being found in the set of cell identifiers pre-associated with the predetermined region; and wherein the location having the set of coordinates is determined to be not within the predetermined region in response to the match not being found in the set of cell identifiers pre-associated with the predetermined region.
 4. The system of claim 3, wherein the coordinates are combined to provide the identifier of the cell according to a predetermined function of the coordinates of the respective location.
 5. The system of claim 3, wherein the coordinates are combined by: generating two integers from the coordinates according to a precision level; and combining the two integers to provide the identifier of the cell.
 6. The system of claim 3, wherein the coordinates are combined by: selecting digits from the coordinates in accordance with a cell resolution level; and combining the digits selected from the coordinates into an integer representing the identifier of the cell.
 7. The system of claim 1, wherein the first mobile devices are identified to have the first status for having received the message during the time period; and the second mobile devices are identified to have the second status for having not received the message during the time period.
 8. The system of claim 7, wherein the subset of the plurality of mobile devices is selected by excluding mobile devices in one or more first clusters that contain only mobile devices having the second status.
 9. The system of claim 8, wherein the subset of the plurality of mobile devices is selected by further excluding mobile devices in one or more second clusters that contain only mobile devices having the first status.
 10. The system of claim 9, wherein in at least a first cluster of the plurality of clusters, devices in the first mobile devices which are also in the first cluster having a population size that is different from a population size of devices in the second mobile devices which are also in the first cluster.
 11. The system of claim 9, wherein a population size of the first mobile devices is different from a population size of the second mobile devices.
 12. The system of claim 11, wherein for each respective cluster that contains at least a portion of the first mobile devices and the second mobile devices, the computing device computes: a first average of responses of devices that are in both the first mobile devices and in the respective cluster; a second average of response of devices that are in both the second devices and in the respective cluster; and a difference between the first average and the second average as a measurement for the respective cluster.
 13. The system of claim 12, wherein the measurement of the location behavior change is computed from a weighted average of measurements computed for clusters that contain the first mobile devices and second mobile devices.
 14. The system of claim 13, wherein the weighted average is based on weights assigned to the clusters that contain the first mobile devices and second mobile devices; and the weights are proportional to sizes of the clusters that contain the first mobile devices and second mobile devices.
 15. The system of claim 13, wherein the sizes of the clusters that contain the first mobile devices and second mobile devices are computed based at least one of: counts of devices that are in both the first mobile devices and the clusters that contain the first mobile devices and second mobile devices; counts of devices that are in both the second mobile devices and the clusters that contain the first mobile devices and second mobile devices; and counts of devices that are in either the first mobile devices or the second mobile devices and in the clusters that contain the first mobile devices and second mobile devices.
 16. The system of claim 7, wherein the subset of the plurality of mobile devices is selected by further selecting from each cluster having a first number of mobile devices having the first status and a second number of mobile devices having the second status; and the first number is equal to the second number.
 17. The system of claim 16, wherein the measurement of the location behavior change is computed by: computing a first average of the responses of first mobile devices; computing a second average of the responses of the second mobile devices; and computing the measurement of the location behavior change from a difference between the first average and the second average.
 18. A method implemented in a computing device, the method comprising: receiving location information of a plurality of mobile devices, wherein each mobile device in the plurality of mobile devices determines locations of the mobile device at respective time instances when the mobile device is at the respective locations; receiving information identifying a time period during which at least one message related to a predetermined region is transmitted to a portion of the mobile devices at least once; identifying a plurality of clusters of the plurality of mobile devices according to a set of attributes, wherein mobile devices in each of plurality of clusters have values of the attributes that are close to each other in an attribute space defined by the set of attributes; selecting a subset of the plurality of mobile devices based at least in part on identifications of the clusters in the attribute space; identifying a plurality of first mobile devices in the subset, wherein each mobile device of the first mobile devices having a first status with respect to the message; identifying a plurality of second mobile devices in the subset, wherein each mobile device of the second mobile devices having a second status with respect to the message; for each respective mobile device in the first mobile devices and the second mobile devices, computing a response of the respective mobile device based on counts of visitation to the predetermined region by the respective mobile device; and calculating a measurement of a location behavior change based on a difference between responses of the first mobile devices and responses of the second mobile devices.
 19. The method of claim 18, wherein the calculating comprising: for each respective cluster that contains at least a portion of the first mobile devices and the second mobile devices, computing a cluster-level measurement of the location behavior change from a difference between: a first average of responses of devices that are both in the plurality of the first mobile devices and in the respective cluster; and a second average of responses of devices that are both in the plurality of the second mobile devices and in the respective cluster; and combining cluster-level measurements of the location behavior change of clusters that have mobile devices in the subset to generate the measurement of the location behavior change.
 20. A non-transitory computer storage medium storing instructions which when executed by a computing device, cause the computing device to perform a method, the method comprising: receiving location information of a plurality of mobile devices, wherein each mobile device in the plurality of mobile devices determines locations of the mobile device at respective time instances when the mobile device is at the respective locations; receiving information identifying a time period during which at least one message related to a predetermined region is transmitted to a portion of the mobile devices at least once; identifying a plurality of clusters of the plurality of mobile devices according to a set of attributes, wherein mobile devices in each of plurality of clusters have values of the attributes that are close to each other in an attribute space defined by the set of attributes; selecting a subset of the plurality of mobile devices based at least in part on identifications of the clusters in the attribute space; identifying a plurality of first mobile devices in the subset, wherein each mobile device of the first mobile devices having a first status with respect to the message; identifying a plurality of second mobile devices in the subset, wherein each mobile device of the second mobile devices having a second status with respect to the message; for each respective mobile device in the first mobile devices and the second mobile devices, computing a response of the respective mobile device based on counts of visitation to the predetermined region by the respective mobile device; and calculating a measurement of a location behavior change based on a difference between responses of the first mobile devices and responses of the second mobile devices. 